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ABSTRACT 



, Context. Kinematical and time variability studies of protostellar jets are fundamental for understand the dynamics and the physics of 

CL(' these objects. Such studies remain very sporadic, since they require long baselines before they can be accomplished. 

Aims. We present for the first time a multi-epoch (20 years baseline) kinematical investigation of HH52, 53, and 54 at optical and 
near-IR wavelengths, along with medium (optical) and high resolution (NIR) spectroscopic analyses, probing the kinematical and 
physical time variability conditions of the gas along the flows. 



o 



^ ■ Methods. By means of multi-epoch and multi-wavelength narrow-band images, we derived proper motions (P.M.s), tangential veloc- 

^ I ' ities, velocity and flux variability of the knots. Radial velocities and physical parameters of the gas were derived from spectroscopy. 

Finally, spatial velocities and inclination of the flows were obtained by combining both imaging and spectroscopy. 
Results. The P.M. analysis reveals three distinct, partially overlapping outflows. Spatial velocities of the knots vary from 50kms"' 
to 120 km s"'. The inclinations of the three flows are 58+3°, 84±2°, and 67+3° (HH52, HH53, and HH54 flows, respectively). In 
20 years, about 60% of the observed knots show some degree of flux variability. Our set of observations apparently indicates accel- 
I eration and deceleration in a variety of knots along the jets. For about 20% of the knots, mostly coincident with working surfaces 

1^-^ , or interacting knots along the flows, a relevant variability in both flux and velocity is observed. We argue that both variabilities are 

. related and that all or part of the kinetic energy lost by the interacting knots is successively radiated. The physical parameters derived 

• ' from the diagnostics are quite homogeneous along and among the three outflows. The analysis indicates the presence of very light 

(TVh ~ 10^ cm '), ionised (X^ ~0.2-0.6), and hot (T^ ~ 14 000-26 000 K) flows, impacting a denser medium. Several knots are deflected, 
, especially in the HH 52 flow. At least for a couple of them (HH 54 G and GO), the deflection originates from the collision of the two. 

For the more massive parts of the flow, the deflection is likely the result of the flow collision with a dense cloud or with clumps. 
' Finally, we discuss the possible driving sources of the flows. 
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1 . Introduction lines in the near-infrared (NIR) (see e. g. Schwartz et al. 119871 

EisloflFel et al. |1994a| Davis et al. |1994| l. At later phases - Class 

Protostellar jets and outflows represent fundamental processes j ^lass II - the sources also become visible in the opti- 

m star formation. These outflows are a direct consequence of ^al, and their jets can then be traced by optical lines as [Sii] 

the accretion inechanism in young stellar objects (YSOs) dmng ^^^^^^ ^^^qJ^^ ^^^^ ^ ^^.^j^^^ ^ Elias[T983l 

their earliest phase (see e. g. Cabrit et al. 1990; Shu et al.|2000l ^^^^^ [T998bt . Observations of the emission lines from the 

Reipurth & Bally 2001), when a traction ot the intalling ma- uiu.j n »»ju -i » j 

. . , ■ , n n . • r„, . . shock-heated gas allow us to study physical parameters m de- 

terial is elected in the form of collimated lets. The interaction . ., , -.u j »u » / »u • ■ e ^■ 

i . , , . „ ,. . tail, such as the density, the temperature, the lonisation fraction, 

between the electa and the circumstellar medium occurs via ra- a u a i u ^ i nnn/i d • * 

^ ^ h k and the abundances (see e. g. Haitigan et al. [1994; Bacciotti & 

Eisloffel 1999; Lavalley-Fouquet et al. 120001 ). Moreover, since 

In the very young and still deeply embedded sources, these ^j^^^,^^ outflows trace ejecta fliat are progressively older with 

flows are mostly observed in molecular (H2) and atomic ([Fe 11]) ^j^j^jj^g f^^m the source, it is also possible to gain important in- 

formation about their driving source properties, reconstructing 

Send offprint requests to: A. Caratti o Garatti e-mail: the YSO mass ejection history. 
alessio@cp . dias . ie 

* Based on observations collected at the European Southern ^ description of the gas dynamics, however, only be- 

Observatory, (Pai'anal and La Silla, Chile) (55.C-0908, 62.I-0136(A), comes available through measurements of the flow velocity. 

63.I-0031(A), 69.C-0269(A), 74.C-0235(A), 77.C-0047(B), 078.C- Spectroscopic radial velocity measurements are necessary, but 

0089(A)) not suflicient for deriving the flow speed, since they deliver only 
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the projected velocities along the line of sight. In order to deter- 
mine the true 3-dimensional kinematics, i. e. the absolute flow 
speed and the inclination of the flow with respect to the plane 
of the sky, proper motion measurements of the knots are also 
necessary. With this information in hand, a number of important 
questions can be addressed straightforwardly. 

i) The sources of the outflows can be unambiguously identi- 
fied and knots can be correctly assigned to their driving sources. 

ii) Using multi-wavelength data, molecular, and atomic, it is 
possible to study the different kinematics of the flow by means 
of tracers that indicate various shock conditions and velocities 
of the gas. 

iii) The internal kinematics of the flows can be studied in de- 
tail. Possibly, acceleration, and deceleration can be derived from 
a multi-epoch analysis of the motions and their origin can be 
investigated in detail. 

iv) Flux variability can be studied and related to the motion 
variability in the flow. 

v) Combining the kinematical results and the gas physical 
parameters, it is possible to infer the mass-loss rate and the en- 
ergy budget along the flow. In turn, we can compare these quan- 
tities with the physical properties of the exciting source to inves- 
tigate how it evolves with time. 

To address these points, we collected a multi-epoch and 
multi-wavelength (optical and NIR) data archive of several 
Herbig Haro (HH) objects and protostellar jets. The database 
includes both imaging and spectroscopy (low and high resolu- 
tion), mainly covering the optical atomic (Ha, [S ii]) and NIR 
molecular emission (H2). For the first time it is possible to study 
the kinematics of the outflows during an elapse of time up to 
20 years at different wavelengths, observing the formation and 
the evolution of the single structures inside the jets. This paper 
represents the first work of our project. 

The structure is as follows. In Sect. 2 we shortly describe the 
investigated region. In Sect. 3 our observations are presented. 
Section 4 reports an overview of our results, including morphol- 
ogy of the flows, their kinematics, flux and velocity variability, 
physical parameters of the gas, and mass flux rates. An investi- 
gation on the possible exciting sources of the flows is presented, 
as well. In Sect. 5 a discussion is proposed in terms of kinemat- 
ics vs. flow variability, flow deflection mechanisms, length, and 
dynamical age of the flows. Our conclusions are summarised in 
Sect. 6. Finally, a detailed description of our results on the single 
knots is given in the Appendices. 



2. The investigated region 

HH52, 53, and 54 are located at the northern edge of the 
Chamaleon II (Cha II) cloud (Schwartz I1 977| l at a distance of 
178+18 pcQ (Whittet et al. flWTl i. As shown in Fig.[T] the three 
HH objects are roughly placed along a straight path, moving 
from SW to NE, with a R A. of -55°. The morphology of the HH 
objects is approximatively the same in the Ha and [S 11] filters. 
HH 52 has a bow shock shape pointing towards HH 54. About 
20" NE, a chain of bright knots (HH53C, A, and B), roughly 
aligned E-W, is observed. Moreover, farther NE (~1'), a clus- 
ter of several bright knots (HH 54) is detected. This group spans 
some tens of arcseconds and elongates towards SSW, delineat- 
ing the so-called HH54 streamer (see Sandell et al. 119871 ). In 



this work, they described the morphology of the HH object^ in 
some details, also for the first time reporting a contour map of the 
molecular hydrogen emission in the NIR. The HHs were spectro- 
scopically observed in the optical by Schwartz & Dopita ( 198C| 
and Graham & Hartigan (1988), who found only blue-shifted 
radial velocities ranging from -40 to -100km s"'. Knee d 19921 ) 
detected two blue-shifted CO outflows, coinciding with HH 52 
and 53, and HH 54 emissions, and only one red-shifted lobe, lo- 
cated approximately between the two. H2 and [Fe 11] imaging to- 
gether with low-resolution spectra of HH 54 were published by 
Gredel ( |1994| l. The morphology in the H2 and [Feii] filters ap- 
pears similar to the optical images. More low-resolution spectra 
in the NIR, with a detailed study of the excitation of the gas 
in HH 54, have been published by Caratti o Garatti (I2006I I and 
Giannini et al. (12006^ From these authors H2 medium-resolution 
spectra of HH54 and radial velocities are also available. These 
data cannot explain, however, the kinematics of the object. An 
optical RM. analysis of the HHs with a seven year baseline was 
attempted by Schwartz et al. ( 119841 ) without conclusive results, 
because of small RM. values and relatively large errors. 

Finally, the exciting sources of the HHs still remain unclear, 
even after a proper motion analysis of the infrared knots (Caratti 
o Garatti et al.^06), the publication of MI PS maps and YSO 
catalogues for this region (Young et al. 120051 Porras et al. 120071 
Alcala et al. lIOOSl l. 

3. Observations, data reduction & analysis 

Our multi-wavelength database is mainly composed of optical 
and NIR observations, collected with the ESO facilities during 
the past twenty years. The relevant information for both imag- 
ing and spectroscopy is reported in Tables[r|and|2l respectively. 
These data were used to derive the kinematics and the physi- 
cal properties of the studied outflows. The advantage of using a 
multi-epoch image dataset is quite straightforward: on one hand, 
we can reduce the errorbars of the RM.s and check for veloc- 
ity variations in the flows; on the other, we can study their flux 
variability and their evolution with time. 

Additionally, this paper makes use of Spitzer Space 
Telescope public data from the "Cores to Disks" (C2d) Legacy 
project (Evans et al. 20031 ), already published and analysed 
(Young et al. 2005 , Porras et al. lIOOTl Alcala et al. .2008j in order 
to infer the exciting sources and map the colder gas component 
of the outflows. 

3.1. Imaging 
3.1.1. Optical 

Our optical narrow-band images (see Table [1]) are centred on 
Ha and [S 11] emissions. A first set of data was obtained at 
the ESO/MPG2.2-m telescope (between 1987 and 1993) with 
three diff^erent detectors: an RCA 30yum CCD chip (1987) (2 
images) and an RCA 15yum CCD chip (1989-1990) (4 images), 
both mounted on the Cassegrain focus, and EFOSC2 (Buzzoni et 
al. 1984 ) (1991-1993) (4 images). A second set of images was re- 
trieved at the New Technology Telescope (1995 and 2006) with 
EMMI (Dekker et al. |1986K 3 images). Because of the relatively 
smafl FoV of the CCDs used between 1987 and 1990, HH52 
and HH53, and HH54 were observed separately with the [S 11] 
(AA 6716, 6731) and with the Ha filters, respectively. Starting 



' Hereafter we refer to this value to compute spatial distances and 
tangential velocities of the knots. 



^ The nomenclature used in this paper mainly refers to Sandell et 

al. iiwn . 
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Date of obs. 


Telescope/ 


Filter 


Resolution 


seeing 


Exp. Time 


Notes 


(d,m,y) 


Instrument 


Band 


("/pixel) 


(") 


(s) 


Objects Observed 


14.05.2006 


ESO-NTT/EMMI 


Ha 


0.166 


0.8 


3600 


HH52,HH53,HH54 


12.05.1995 


ESO-NTT/EMMI 


Ha 


0.266 


1.1 


3600 


HH52,HH53,HH54 


31.03.1993 


ESO/MPG2.2-m/EFOSC2 


Ha 


0.35 


0.9 


4500 


HH52,HH53,HH54 


08.05.1992 


ESO/MPG2.2-m/EFOSC2 


Ha 


0.35 


1.4 


5400 


HH52,HH53,HH54 


15.01. 1990 


ESO/MPG2.2-m/CCD 


Ha 


0.35 


1.0 


6000 


HH54 


03.02.1989 


ESO/MPG2.2-m/CCD 


Ha 


0.35 


1.0 


4200 


HH54 


04.03.1987 


ESO/MPG2.2-m/CCD 


Ha 


0.35 


1.2 


3600 


HH54 


14.05.2006 


ESO-NTT/EMMI 


[Sii] 


0.166 


0.7 


3600 


HH52,HH53,HH54 


31.03.1993 


ESO/MPG2.2-m/EFOSC2 


[Sii] 


0.35 


1.0 


4200 


HH52,HH53,HH54 


09.02.1991 


ESO/MPG2.2-m/EFOSC2 


[Sii] 


0.35 


1.1 


1800 


HH52,HH53,HH54 


14.01. 1990 


ESO/MPG2.2-m/CCD 


[Sii] 


0.35 


0.8 


4800 


HH52,HH53 


01.02.1989 


ESO/MPG2.2-m/CCD 


[Sii] 


0.35 


1.2 


3600 


HH52,HH53 


28.02.1987 


ESO/MPG2.2-m/CCD 


[Sii] 


0.35 


1.5 


3600 


HH52,HH53 


01.01.2005 


ESO-VLT/ISAAC 


H2 


0.148 


0.5 


40 


HH54 


02.07.2002 


ESO-VLT/ISAAC 


H2 


0.148 


0.5 


60 


HH54 


06.06.1999 


ESO-NTT/SofI 


Hz 


0.29 


1.0 


300 


HH52,HH53,HH54 


31.03.1999 


ESO-NTT/SofI 


H2 


0.29 


0.9 


600 


HH52,HH53,HH54 


15.03.1995 


ESO/MPG2.2-m/IRAC2b 


H2 


0.51 


1.1 


500 


HH52,HH53,HH54 


10.04.1993 


ESO/MPG2.2-m/IRAC2b 


H2 


0.51 


1.7 


200 


HH54 


03.04.1993 


ESO/MPG2.2-m/IRAC2b 


H2 


0.51 


1.7 


200 


HH52,HH53 


01.09.2007 


ESO-NTT/SofI 


H2 


0.29 


1.0 


1080 


IRAS 12500-7658 


01.09.2007 


ESO-NTT/SofI 


K, 


0.29 


1.0 


270 


IRAS 12500-7658 



Table 2. Journal of observations - Spectroscopy 



Telescope/ 


Date of obs. 


Wavelength 




PA. 


K 


Encompassed Knots 


Instrument 


(d,m,y) 


(A) 


(s) 


(°) 






ESO-NTT/ 


14.05.2006 


6200-7000 


2700 


52 


3000 


HH52: A3-A4,AI,D4; HH53: C-C1,E1; 


EMMI 












HH54: G-G0,G1-G3,C2,C3 


ESO-VLT/ISAAC 


01.01.2005 


16400 


3600 


29.5 


10000 


HH54: Z,A3,A1,BJ-J1,H3,H2 


ESO-VLT/ISAAC 


01.01.2005 


16400 


3600 


160.5 


10000 


HH54:C3,C1,H3,J-J1,E,K 


ESO-2.2-m/ 


20.02.1987 


6200-6800 


3600 


89.5 


2000 


HH53: B,A,C,F2; HH54: X4C,X1 D 


Boiler & Chivens (B&C) 


21.02.1987 


6200-6800 


3600 


89.5 


2000 


HH53: B,A,C; HH54: X3,X1 A 




25.02.1987 


6200-6800 


3600 


75 


2000 


HH52: B,A,D3-D2; HH54: X4B-A 



from 1991, all the three objects were observed in both filters 
within a single frame. The exposure times range from 1800 to 
6000 s, resulting in deep images with a high signal-to-noise ratio 
(S/N). The raw data were reduced by using standard procedures 
in //M/Q, subtracting bias, flat-fielding and removing bad pixels. 
Cosmic rays were removed by using the L.A. Cosmic algorithm 
(van Dokkum |20011 l. which was extremely efficient in cleaning 
cosmic-ray crowded images. The measured seeing on our im- 
ages ranges between (X'T and I'.'S (Table[T]), with the best values 
obtained in the most recent ones. 

Only images from 2006 have their own photometric standard 
star, observed with both Ho; and [S ii] filters. The selected stan- 
dard is LTT 3218 (Hamuy et al. |1994| l. for which a flux calibrated 
spectrum also exists. At the required wavelengths, the flux of the 



^ IRAF (Image Reduction and Analysis Facility) is distributed by 
the National Optical Astronomy Observatories, which are operated 
by AURA, Inc., cooperative agreement with the National Science 
Foundation. 



star was obtained by convolving its flux density with the profiles 
of the ESO filters. 

To flux-calibrate the earlier epoch images, 30 bright field 
stars were selected and photometry in both filters performed on 
them. For each epoch in each filter, the star fluxes were cross- 
checked, and those showing variations less than 3% were chosen 
as secondary standard stars to calibrate the fields. 

3.1.2. NIR 

Our narrow-band H2 (2.12yum) images were collected during 
different runs between 1993 and 2005 , at the ESO/MPG2.2-m 
telescope with IRAC2b (Lidman et al. |1998al l (3 images), at the 
ESO-NTT with Sofl (Moorwood et al. 199831 (2 images), and at 
ESO-VLT with ISAAC (Moorwood et al. 1998b) (2 images) (see 
Table [TJ. Some data on HH54 have been presented in previous 
papers (Giannini et al. 2006; Caratti o Garatti et al. 2006), where 
P.M.s were derived comparing the Sofl image of June 1999 with 
the ISAAC one of January 2005. Here, we enlarge the baseline 
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I I I I I I I I I I I I 1 I I I I I 1 I I 1 I I , I i_ 

a).0 40.0 30.0 20.0 10.0 12:55:00.0 

R A. (j200(}) 



Fig. 1. [S ii] EMMI 2006 image of the region indicates the locations of the HH objects, along with EMMI, B&C, and ISAAC slit 
positions. In the upper-right inset [S ii] deep image (from 1993) we show the newly detected HH 52 streamer. A rectangular box (in 
the centre of the image) shows its exact location. Additionally, left and right boxes approximatively indicate the fields displayed in 
Figuresglll] and|6] 



of our observations including IRAC2b data from 1993 and 1995, 
and an ISAAC image from 2002 (covering only HH 54) that, on 
one hand, allow us to also derive H2 P.M.s for HH 52 and HH 53, 
and, on the other, add more data points to HH 54, reducing the 
uncertainties on P.M.s and P.A.s of the knots. The first set of 
IRAC2b data (1993) covers HH52 and HH53, and HH54, sep- 
arately, exhibiting a poor seeing of 1'.'7. The second (1995) con- 
sists of a mosaic of the region (~7'x7') and has a seeing of I'.'l. 
In addition, more H2 SofI data from the ESO science archive fa- 
cilitjQ (taken less then 3 months apart from those of June 1999, 
see Table [TJ were retrieved and coadded to create a deeper SofI 
H2 map of the region covering ~irxir. 

Finally, additional H2 and K, images around IRAS 12500 - 
7658, one of the possible driving sources of HH54 (see also 
Sects. [3T3] and \4^ . were taken in Jan. 2007 to detect a pos- 
sible H2 jet from that source. 

As for the optical, all the raw data were reduced by using 
IRAF packages applying standard procedures for sky subtrac- 
tion, dome flat-fielding, bad pixel and cosmic ray removal and 
image mosaicing. Flux calibration by means of a photometric 
standard star was only possible for the June 1999 SofI image. 
As for the optical, the remaining H2 images were calibrated by 
selecting a few field stars that showed flux variations smaller 
than 10%. 



http://archive.eso.org/ 



3.1.3. Spitzerdata 

Reduced Spitzer data have been obtained by the public catalogue 
c2d (third delivery). They consist of large mosaics of the Cha II 
cloud observed with the Infrared Array Camera (IRAC) in four 
channels (at 3.6, 4.5, 5.8, and 8.0 pm) (~1.04deg2) and with 
MIPS (at 24, 70 yum) (-1.5 deg^). In Fig. |2] we report part of 
the Spitzer-MIPS map (at 24yum), indicating the location of the 
studied HH objects and the possible driving sources. 

3.1.4. Knot identification and flux analysis 

To identify the knots along the flows, we used the best image 
in each filter (i. e. with highest S/N and best resolution). This 
resulted to be the last epoch image (i. e. EMMI 2006 image in 
Hff and [S 11], SofI 1999 image for HH 52 and HH 53 in H2, and 
ISAAC 2005 image for HH54 in H2). These were also used as 
reference images for the P.M. analysis (see Sect. I3.1.5I I. Each 
structure was then identified and labelled (in each filter) within 
a 5 cr contour around the local maximum. These contours were 
also used to measure the total flux of the knots in each epoch and 
filter, and to study the flux variability as well. 

3.1.5. P.M. analysis 

Particular attention was paid to align the images. In fact, an 
accurate determination of the proper motions requires that the 
multi-epoch images have to be scaled and then registered to a 
subpixel accuracy in each filter. Our observed region contains 
several field stars: therefore, the image registration was easy to 
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56:00.0 54:00.0 52:00.0 12:50:00.0 48:00.0 46:00.0 44:00.0 

R.A. (J2000) 

Fig. 2. Spitzer-MIPS 24fj.m map of the investigated Cha II region. The candidate exciting sources and the HH objects are reported. 
A rectangular box indicates the region displayed in Fig.[T] 



implement, and the distortions introduced by the different optical 
systems were corrected by means of geomap (with a polynomial 
fit of 3rd order) and geotran routines in IRAF. Because of the 
large baseline of the observations, few stars showed a significant 
proper motion and were excluded from the sample. Each image 
in each filter was aligned to a common (scaled) reference frame 
(see Sect. l3.1.4l i. The resulting error is given by the residuals of 
the fit. Statistical errors in the optical were typically about +0. 1- 
0.2 pixels (i. e. ~0'.'03-0'.'07), depending on the seeing and on 
the resolution of the image. In the NIR, due to the reduced FoV, 
ISAAC, and IRAC2b (1993) images have an exiguous number 
of stars (<20) that diminishes the accuracy in the alignment, not 
allowing the entire area to be sampled homogeneously. As a con- 
sequence, here the errors were up to ±0.3-0.4 pixels. 

Knot shifts were determined between image pairs (i. e. last 
epoch and previous epoch images) using a cross-correlation 
method. A rectangle was defined round each knot, enclosing its 
5 cr contour. The 1 st-epoch image was then shifted with subpixel 
accuracy (0. 1 pix) with respect to the last epoch image and then 
multiplied by it. For each shift (x,y), the total integrated flux (f) 
in the rectangle around each knot was measured in the product 
image. As the final shift for each knot, we used the position of 
the maximum of f(x,y), determined via a Gaussian fit. For each 
structure, an initial guess of the shift vector was estimated vi- 
sually by blinking the image pairs, defining a maximum shift 
for each rectangle. This served as input into our routine for the 
cross-correlation. For several features, the analysis routine con- 
verged rapidly to a solution, independent of the initial guess. A 
more accurate estimate (within a few pixels) of the size and di- 
rection of the proper motion was needed to avoid the code to 
correlate spurious pairs in case of diffuse knots located close to 
brighter knots or stars. To prove the code accuracy, we tested it 
on both nebulous and stellar objects by shifting the images and 
recovering the offsets with the code. 



2.5 

2.0 

X 1.5 
'q. 

Z 1.0 

"5 0.5 

X 

0.0 



HH54Y2 




1995 2000 
Epoch (yr) 




1995 2000 
Epoch (yr) 



Fig. 3. Example of linear and quadratic fits performed on the 
dataset to derive proper and accelerated motions. Here the fits 
on HH 54 Y2 Ha data are reported. 



To quantify the systematic errors for the shift measurements, 
we varied the size and shape of the rectangle that defines the 
object. The resulting range of values gives a systematic error, 
which depends on the S/N, shape, and time variability of the 
knot. This error is comparable to or even larger than the align- 
ment error Thus, each single shift error was derived combining 
both the alignment error of the two epoch images and the un- 
certainty in the centring routine. Usually, for each pair, the final 
error was around 10% of the P.M. value. However, for slow knots 
with a low S/N, the error was up to 100%. 
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The final P.M. value and error for each knot in each filter 
were obtained from a weighted least square fit of the shifts, 
fitting the motion in X (R. A.) and Y (Dec.) simultaneously 
(see Fig. |3]l, using the equations of both linear and accelerated 
motions: [x(t),y(t)] = p^y + v^yt + {jay-yt^). Both linear and 
quadratic fits (for knots with a number of epochs > 3) were per- 
formed to detect any possible acceleration or deceleration of the 
knot, as well. The resulting errors on proper and accelerated mo- 
tions were derived from the fits. An example of our fits is shown 
in Fig. |3] where an accelerated motion is detected. The upper 
and lower panels report the X and Y shifts with errors for the 
Ha filter of HH 54 Y2, respectively. Straight and dashed lines 
indicate linear and quadratic fits. In this particular case both X 
and Y shifts are better modelled by a quadratic fit. 

3.2. Spectroscopy 

3.2.1. Optical 

Four optical medium-resolution spectra CR ~2000-3000) were 
collected with a Boiler & Chivens (B&C) spectrograph at the 
ESO/MPG2.2-m telescope (1987) and with EMMI at the NTT 
(2006). They range from 6200 to 7000 A and cover the spectral 
region containing bright Unes from [O i] (6300 A and 6363 A), 
[Nn] (6548 A and 6583 A), Ho- (6563 A), and [Sn] (67 16 A and 
6731 A), widely used for the diagnostics of the shocked gas (see 
e.g. Bacciotti & Eisloffel [19991 hereafter BE99). The observed 
targets and the slit P. A. are listed in Table |2l together with the 
date of the observations and the total integration time. B&C slit 1 
and 2 are positioned on the three HH53 bright knots (A, B and 
C) (see also Fig.[T]), and intersect the HH 54 streamer They have 
the same PA. (89°5) and almost the same positions, with slit 1 
shifted slightly northward and also encompassing knots F2 in 
HH53 and X4C and XI D in HH54, while slit 2 covers knots 
X3 and XI A of HH54 streamer (see Fig.|4]and AppendixlAlfor a 
detailed description of the flow morphology). B&C slit 3 (PA. - 
75°) is placed on HH 52 A, B and D2-D3, and also includes knots 
X4 A and B in HH54 streamer. Finally, EMMI sUt (PA. = 52°) 
covers HH52B, A, D3-D2, HH53C and CI, HH54G, and C. 
The total integration time for B&C and EMMI spectra are 3600 s 
and 2700 s, respectively. We used the same procedure as for the 
optical imaging to reduce the data. A first rough wavelength cal- 
ibration was retrieved from helium and argon lamps. A second, 
more accurate calibration was done using several atmospheric 
OH lines (Osterbrock et al. 1996) present in each frame with a 
final accuracy of 4 km s"' in the EMMI spectrum. LTT 7379 was 
the standard star observed to flux calibrate the EMMI spectrum. 
The B&C data were calibrated using our narrow-band images. 

3.2.2. NIR 

Two medium-resolution spectra around the 1 .64 jum [Fe ii] line 
i'^vac = 1. 64399 //m; Johansson 1978) were obtained with the 
ISAAC IR camera, with a 0'.'3 slit, corresponding to a nominal 
resolving power of ~ 10 000. Targets and P.A.s of the slits are 
reported in Table |2l while the slit positions are shown in Fig.[T] 
(ISAAC slit 1 and 2). The two slits encompass only HH54. The 
first slit (PA. = 29^5) encompasses knots A3, Al, B3, Bl, par- 
tially J and Jl, H3, and H2. The second sHt (PA. = 160?5) is 
centred on knots Kl, K, and E, but partially encompasses knots 
J, Jl,H3,Cl,and C3. 

To perform our spectroscopic measurements, we adopted 
the usual ABB'A' configuration, for a total integration time of 
3600 s. Each observation was flat-fielded, sky-subtracted, and 



corrected for the curvature derived by longslit spectroscopy, 
while atmospheric features were removed by dividing the spec- 
tra by a telluric standard star (08 spectral type), normalised 
to the blackbody function at the stellar temperature and cor- 
rected for hydrogen recombination absorption features intrinsic 
to the star The wavelength calibration was performed using at- 
mospheric OH lines (Rousselot et al. 2000) with an accuracy 
of 3 km s ' . The instrumental profile in the dispersion direction, 
measured from Gaussian fits to the sky OH emission lines, was 
38.0kms-i. 



4. Results 

4.1. Morphology 

An overview of the HH 52, 53, and 54 reg ions is given in Fig.[T] 
Additionally, in Figures |4l |A. U and lA.jfl we present Ha, [Sn], 
and H2 images of HH 52, HH 53, and HH 54 with labels indicat- 
ing the position of the knots, including the newly detected ones. 
Thei r coordinates are reported in Columns 2 and 3 of Tab. IB.U 
IB.4fI Note that some emissions are not visible in all the three 
filters. The atomic and molecular emissions only partially over- 
lap, sometimes showing a different spatial distribution of the gas 
in each filter that indicates various excitational conditions along 
the flows. 

A first glance to Fig. [T|reveals at least two flows converging 
towards HH 54 brightest region, which shows a chaotic struc- 
ture (see also Fig. |4|i. A first, well-known outflow, delineated by 
the HH54 streamer, is oriented NNE with a position angle of 
~22° (see also Fig. HJi. A second outflow follows an NE direc- 
tion with a position angle of 50°-55°, grouping HH52, 53, and 
54. In our images we detect a new streamer (hereafter HH 52 
streamer, see upper-right panel in Fig. [1]), which connects the 
HH 52 bow shock to the main body of HH 54 and it is super- 
imposed on HH53. The structure is barely visible in H2, but 
well delineated in our deepest [S 11] and Ha images. It is de- 
tected in the Spitzer images as well. It appears as continuous 
slightly curved emission between HH54 and 53 knots, while it 
is fragmented in single knots between HH53 and 52. Both the 
narrow-band images and the EMMI spectrum seem to indicate 
that the HH 52 streamer only presents line emission. 

We do not detect any other additional emission in the 
large H2 Soft map of the region down to the 3cr limit 
of ~3xlO"'^ erg s"' cm"^ arcsec"^. Neither do we detect 
any jet or molecular emission in the 2.12//m image around 
IRAS 12500 - 7658 (a possible driving source of HH54 
- see also Fig. |2] and Sect. |4.9l l down to a 3cr limit of 
~1.5xlO"'^ erg s"' cm"^ arcsec"^. Moreover, we do not ob- 
serve any further emission in the Spitzer large field images, with 
the exception of a faint nebulosity in IRAC band 2 (4.5 ^m) lo- 
cated south ~2'.8 SSE of HH52 (at about a = 12''54'"50' and 
S =-77°00'37'.'2) and extending approximatively ~\' towards 
IRAS 12496-7650, located ~8' farther SSE. It is not clear, how- 
ever, whether this emission is correlated with the outflows. 

Finally, in the Appendix |A] for each object, a detailed de- 
scription of the morphology is reported. 



' Figs. |A.l1 and |A.2| are presented in Appendix [A] as online material 
only, available in electronic form at http ://ww w. edpscience s . org [ 
* Tables ls! 1 IB .4l are presented in AppendixiBlas online material only. 
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Fig. 4. Ua images (EMMI 2006) of HH54 (left panel), HH52 and HH53 (right panel). The labels indicate the position of the 
knots, including the newly detected ones. The contour levels are 3, 10, 20, 30, 40, 50, 60, 70, 80, 90, lOOx the standard deviation to 
the mean background (cr ~ 4 x 10"'^ erg s"' cm"^ arcsec"^). 



4.2. Proper motions 

In this section we report the general results of our proper motion 
analysis. A detailed description for each HH object is given sep- 
arately in Appendix |B] During the time interval of our observa- 
tions, some objects showed high degrees of variability, consider- 
ably changing in velocity, direction, shape, and flux. Sometimes 
it was possible to observe the formation or the vanishing of struc- 
tures along the flows. Therefore, we leave to Sects. 143114.41 and 
to Appendix |C] the detailed report of these phenomena. 

InFigs.lSlil andlBllTl we show the flow charts of the two 
regions for the three filters (Ha, [S ii], and H2, respectively). The 
P.M.s and their error bars are indicated by arrows and ellipses, 
respectively. The results of our P.M. measurements (from the lin- 
ear fits) are given as online material in Appendix IB] Table IBTI 
(HH52 and 53) and Table 12] (HH 54) report the values of the 
optical lines, and Table |B3] (HH 52 and 53) and |B4l (HH 54) 
those of the H2 line. Each table lists the knot ID, the coordinates, 
the measured proper motion from the linear fit, the derived tan- 
gential velocity (vtan) at a distance of 178 pc, and the filter used. 
Where P.M.s could not be measured only the knot ID, the coordi- 
nates, and the filter (in which the knot was detected) are reported. 
Those data which are better modelled by a quadratic fit are la- 



belled with a n 'a', and the results of the quadratic fits are listed 
in Table|cl|(see Sec. 1431). 

With few exceptions, the derived P.M.s have similar values in 
the optical filters, with values ranging between 0.007 and 0.097 
" yr"', corresponding to a Vtan between 6 and 82kms"'. On aver- 
age, the H2 P.M.s have lower values, ranging between 0.007 and 
0.041 " yr"' (or 6 and 35 km s"'), and the uncertainties are usu- 
ally greater, especially for the HH52 and 53 regions, due to the 
smaller number of epochs, poor seeing, lower resolution of the 
used instruments, and less accuracy in the alignment (see also 
Sect, mill- 

A quick inspection of HH52 and 53 (right panel of Figs. |5] 
and |6]l reveals the presence of two diff'erent outflows. A first 
flow (hereafter HH 52 flow) is composed of the majority of the 
knots (those of HH52 and several of HH53, i.e. from CI to 
G). It moves NE with a position angle of ~55° (towards HH 54), 
roughly following the path delineated by the HH 52 streamer and 
partially superposing on a second flow (hereafter HH53 flow). 
This one, made of HH53 A, B, Bl, and C, has approximatively 
an E-W direction and does not seem to be related to the HH 52 
flow. The features of the HH52 flow are moving at a Vtan be- 
tween ~30 and 60kms where the different velocities do not 
seem to be associated with the position of the objects along the 
flow, but rather they could originate from local interactions in- 



^ Fig. lB.ll is presented in Appendix IbI as online material only. 



** Table lC. fl is presented in Appendix|C]as online material only. 
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Fig. 5. Flow charts of HH54 (left panel), HH52 and 53 (right panel) in the Ha fiher. Proper motions and their error bars are 
indicated by arrows and ellipses, respectively. 



side the fluid. On the contrary, HH 53 A, B and C, in the HH 53 
flow, exhibit smaller P.M.s (0.01-0.02" yr^' or 9-18 km s"') with 
P.A.s around 270°, indicating that they are probably part of a 
dififerent outflow. 

The left panels of Figs. |5]and|6]reveal the complex kinemat- 
ics of HH 54. The trajectories of the knots indicate that at least 
two distinct outflows are present in this region and they appear 
to be superimposed on the main body of HH54. A first flow 
is generated by the HH 54 streamer {hereafter HH 54 flow) and 
is moving NNE (P.A.~20°). It includes knot groups X, Y and 
Z (HH 54 streamer) and likely groups E, K and Q (in the main 
body), with the Vtan that decreases from ~60 to ~10kms"' mov- 
ing from the stream to the main body. A second outflow with 
an NE direction (P.A.~55°) is the extension of the HH52 flow 
and contains groups G, C, and, possibly. A, which appears de- 
flected (P.A.~80°, Vtan ~30kms"'). Indeed, for several knots the 
affiliation with one of the two flows cannot be certain, and the 
presence of other flows, even if unlikely, cannot be ruled out (see 
also Appendix IC] for a detailed discussion). This is particularly 
evident in the central part of the HH object, where the overlap of 
the flows and the low P.M. values make it difficult to distinguish 
between real and apparent motions. 

4.3. Accelerated motions 

As anticipated in the previous section. Table IC.ll reports those 
datasets that are better modelled by a quadratic fit. In particular, 
the table lists, the knot ID; the goodness of the chi-square for the 
quadratic fit (Paix^/v), i.e. the 1-tail probability value associ- 



ated with the provided chi-square value and degrees of freedom 
v); the proper and accelerated motions derived from the fit; the 
P.A.s of the P.M. and of the acceleration vector; the filter and 
the notes. The data have been selected by following criteria: a) 
P(X^/v) of the quadratic fit is greater than that of the linear fit; 
b) Pi(x^/v)>0.05; c) the measured accelerated motion is >2cr. 
For datasets with only three measurements, conditions a and b 
are always verified with P^ix^ /v)-l. In this case we report the 
result of the quadratic fit only to support or reject the evidence 
of an accelerated motion measured in other filters. It is worth 
to note that same or opposite directions of proper and acceler- 
ated motion (i e. PA. values of v and a, respectively) indicate 
acceleration or deceleration of the knot. As a result of this se- 
lection. Table ICTI lists 33 (out of 108) knots that seem to show 
some velocity variability in at least one filter. However, not all of 
these results seem to be genuine. For example, a few knots (i. e. 
HH 52 A5, 53 G, 54 G2, 54 H3, see Table|OJ have inconsistent 
PA. (a) values in different filters, whereas we would expect sim- 
ilar dynamics for the same knot even with different filters (i. e. 
it is unlikely that the same knot is accelerated in a filter and de- 
celerated in the other). Moreover the P.A.(v) of some other knots 
appears to be inconsistent with respect to the PA. measured from 
the linear analysis, and with respect to the PA. that is expected 
from the geometry of the flow. We labelled all these knots as 'no' 
in the notes of Table lC.ll Indeed for these knots the quadratic fits 
do not represent a reliable physical solution. After this further se- 
lection, only 17 knots are left. They indeed appear to be genuine 
accelerated motions. Among them, 12 objects accelerate and 5 
decelerate. For most of them, both modules and P.A.s of proper 



A.Caratti o Garatti et al.: Time- variability in the Cha II outflows 



9 




and accelerated motions are identical (inside the errors) for the 
two optical filters. 

4.4. Flux, direction, and veiocity variabiiity 

In this section we report some examples of the observed time 
variability in flux, velocity, and direction. In Appendix|C]we pro- 
vide details on those most significant. 

On average, the analysed knot fluxes are variable with time. 
About 60% of the analysed knots (62 out of 108) showed a vari- 
ability above the 3 cr error bars. However, such variations often 
(~60%) do not have a particular trend, but rather appear as tem- 
porary fluctuations around a mean value. For these objects, the 
P.M.s, inside the error bar, are usually constant with time. On 
the other hand, in 20 years about 20% of the knots (i. e. 19 out of 
108) showed velocity variability, together with flux variations. 

The structures that show the highest variability in morphol- 
ogy, flux, and velocity can be mostly identified with (1) working 
surfaces or (2) interacting knots, where the faster fluid particles 
ejected at later times eventually catch up with the slower flow 
ejected earlier from the source. Usually, the fast gas from be- 
hind is decelerated, while the slow gas in front is accelerated. 
This phenomenon is mostly accompanied by flux variability. In 
particular, we identify four regions that clearly illustrate such a 
phenomenon: HH 52 B and C, HH 54 G, and HH 54 H. 

For the analysed knots, the largest flux variations have been 
usually detected in the Ha and H2 emissions, while the [S 11] usu- 
ally presents a lower degree of flux variability. Also the timing 
and the duration of such a variability seem to be different in the 
three filters. In particular, the optical component seems to excite 



first and then raises the molecular emission, which usually seems 
to exhibit a longer cooling time (> 10 years). This is, however, a 
very qualitative analysis, since we do not have H2 measurements 
from the beginning of our time series. Moreover, such a trend can 
be modified by a subsequent shock front. All this can be seen in 
Fig.|2l where three examples of the observed knot flux variabil- 
ity are shown. In addition more examples of flux variability are 
shown in Appendix ICl The measured flux in the dififerent filters 
([S 11] - circles. Ha - triangles, and H2 - squares) is plotted as a 
function of time. The top panel of Fig. [T] shows the appearance 
in 1989 of knot HH52B5, which rapidly increased its bright- 
ness and then faded at optical wavelengths, while the H2 emis- 
sion seems to rise later. Central and bottom panels of the figure 
illustrate two more examples of variability in HH52. HH52C1 
has a trend similar to B5, but the collision with a fast overtak- 
ing shock between 1995 and 2006 (see also Appendix O newly 
increases the Ha and [S 11] luminosity. Finally, knot HH 52 C2 
(bottom panel) represents a quite rare case of [S 11] variability, 
where the knot appears visible in Ha only in 2006. 

In Fig.[8]an example of P.M. and flux variability is indicated. 
The bottom panel gives the P.M. values for the optical filters 
(no molecular emission is detected) of HH 52 B 1 plotted against 
time. As described in Sect. 13.1.51 The P.M. value decreases in 
both the [S 11] and Ha filters (see Table IC.ll column 4), with a 
deceleration of 0.022+0.003 and 0.02+0.01 " yr"-, respectively. 
We also note that the PA.s of the linear motion and of the ac- 
celerated motion (Table IC.ll columns 5 and 6) are identical in 
both species (within the errors). In addition the opposite direc- 
tion of the accelerated motion indicates deceleration of the knot 
for both species. This is also visible in ourlinearfitof Fig.[8](bot- 
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1990 1995 2000 2005 

Epoch (year) 

Fig. 7. Measured fluxes (uncorrected for the extinction) in knots 
HH52B5 (top panel), CI (central panel), and C2 (bottom 
panel) in [S ii] (circles). Ha (triangles), and H2 (squares) filters. 



torn panel), where we consider both species deriving an average 
deceleration value of 0.017+0.005" yr^^ or 0.022+0.002" yr^^ 
without considering the 1990 [S 11] data point, apparently discor- 
dant. Also the flux of the knot varies with time, with a delay of 
about 3 years between the [S 11] and Ha lines (as shown in Fig.[8j 
top panel). 

In Fig.|9] as an example of the observed variability in veloc- 
ity, direction, and flux, the behaviour of optical knots HH 54 GO 
and G is presented. As for the previous knot, our analysis of the 
accelerated motions of these objects can be found in Table ICTl 
These knots are part of the HH52 flow (P.A:~55°) in the NE 
region of the main body of HH54 (see also Fig. [H Fig. IA.2I 
and Appendices A and B), where the fast flow is suddenly de- 
celerated as it collides with the leading material, which is par- 
tially deflected sideways (P.A.~80°). GO is part of such a fast 
flow, that shocks at the end a slow moving component (knot G), 
changing its direction. Bottom panel of the upper figure indi- 
cates the measured P.M.s of GO (circles) and G (triangles) at 
different epochs in both optical filters (open and filled marks 
are Ha and [S 11] data, respectively). The top panel (upper fig- 
ure) shows the variations of P.A.s with time for both knots in 
both filters. Finally, the lower figure reports the flux variabil- 
ity of the knots. The P.M. in GO decreases from -0.07 to -0.04 
"yr"'. (i.e. from a Vtan -65kms"' to a Vtan -35kms"') as it 
interacts with knot G ahead. At the same time, its PA. also ex- 
hibits an abrupt change from ~55° to ~39°. Accordingly, P.M. 
of knot G increases from ~0.03 to -0.06 " yr ' . (i. e. from -25 
to -53kms ' in Vtan), while its P.A. turns from -65° to -80°. 
Indeed such a mechanism could at least partially explain the dif- 
ferent P.A.s observed along the flows. As for the other knots, 
we observe a flux variation in both knots as well (Fig. |9l lower 
figure). 



1990 1995 2000 2005 

Epoch (year) 

Fig. 8. Vai-iability in HH 52 B 1 . Bottom panel shows the P.M.s 
as a function of time in both [S n] (circles), and Ha (triangles). 
The continuous and dashed lines are the best fits of the data 
points with and without 1990 [S 11] point, respectively. The slope 
of the fit gives the deceleration of the knot. In the top panel the 
knot flux (in both filters) (uncorrected for the extinction) as a 
function of time is indicated. 

4.5. Radial velocities 

In Tables IDTIID.^R the results from our spectroscopical analysis 
are reported. From Table lD.ll to lD.3l we list, for each object and 
optical instrument (EMMI and B&C), the detected lines, the vac- 
uum wavelength, radial velocities (vrad), and calibrated fluxes. 
Finally, Table ID . 4l indicates the [Fen] (1.64yum) radial velocities 
along the ISAAC slits. No flux calibration is available for the 
NIR spectra. Each Vrad has been computed with respect to the 
local standard of rest (LSR), assuming a cloud speed of 2 km s"' 
(see Knee 1992 ). When detected, two different velocity compo- 
nents are reported. 

At variance with B&C, the EMMI data allowed us to study 
single substructures inside the HH objects, thanks to the higher 
spatial and spectral resolution, and often two distinct velocity 
components were detected. Moreover, errors in the EMMI spec- 
tra are considerably smaller than in the B&C (4 and 20km s ', 
respectively). However, for those regions encompassed by slits 
of both instruments, there is excellent agreement among the val- 
ues measured with the two instruments. It is also worth noting 
that all the measured velocities of the knots in this region are 
blueshifted, as mentioned by Graham & Hartigan ( I1988I I. With 
a few exceptions, our results are in excellent agreement with 
theirs. Due to reduced resolution, however, both their spectra and 
our B&C spectra only show one velocity component. 

Our measured radial velocities range between -20 and -110 
km s"' with variations similar to those detected in the P.M.s. The 
correlation between radial velocities and P.M.s implies that knot 
velocity variations are real. With a few exceptions, radial veloc- 
ities derived from the detected (atomic) lines have similar val- 

' Tables ID.UD.4I are presented in Appendix |D] as online material 
only. 
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Fig. 9. Upper figure: changes of proper motions (bottom 
panel) and directions (top panel) in HH54G (triangles) and 
GO (circles) with time. Open and filled marks indicate mea- 
surements from Ha and [S ii], respectively. Lower figure: mea- 
sured (uncorrected for the extinction) fluxes in HH 54 G (bottom 
panel) and GO (top panel) in [S ii] (circles), Ha (triangles), and 
H2 (squares) filters 



ues. The highest values (around -llOkms ') are observed in the 
three knots of the HH 53 flow (A, B, C) (see Table \D2^ and 
along the two streamers (around -100km s"'). A change in the 
radial velocity (from ~100kms"' up to ~40kms"') is observed 
ahead of the HH 52 streamer in groups HH 54 G and C, where 
the fast flow impacts slow-moving material. 

The line-profile analysis of the two velocity components in- 
dicates that they originate from single bow-shock structures (as 
in HH52A), from two diff^erent flows (as in HH53C and CI), 



or from a flow at different speeds (e. g. in HH54G) (see details 
in Appendix iDli. 

Figures [To] and [TT] show the position velocity (PV) diagrams 
of the two ISAAC spectra. Since we do not have any [Fe 11] im- 
age, we use an Ha image (EMMI 2006) to indicate the position 
of the slits, enlarged to match the spatial resolution of ISAAC. 
The accordance between the Ha and [Fe 11] emissions is quite 
satisfactory. The observed velocity values range from about - 
10 to -llOkms The velocity structure of the region appears 
complex, and again, often two velocity components are detected. 
At variance with the EMMI spectra,it is more difficult here to 
tell whether these components originate from single bow shock 
structures or from the overlapping of two different flows (i. e. 
HH52 and 54). 

Sometimes the two velocities are well-separated, as in B3 or 
close to B (Fig. [Toll, or as in CI and C3 (Fig.[TTI). For these ob- 
jects, the superposition of the two diff'erent flows is very likely. 
The fast component would then originate from the HH 54 flow, 
which shows in Z a v,-ad of — 100 kms"' (see Fig.[TOli. In other 
cases, such as in Al, B, H2, and H3 along slit 1 (see Fig.[TOli. the 
two velocity components are not well-separated and the line pro- 
files resemble those of a bow shock. For HH 54 B, however, the 
highest intensity in the spectral profile is detected in the lower 
velocity component, as expected in a cloudlet geometry (see, 
e. g., Schwartz il9781 l and not in the higher velocity component, 
as observed in a bullet model (e. g. Hartigan et al. il987. Davis et 
al. l2001l l. In the first case the observed P.M.s should be close to 
zero, and, indeed, in knot B we observe a Vtan close to zero. 



4.6. Inclination and spatial velocity 

By combining tangential and radial velocities it is possible to in- 
fer an upper limit of the inclination of the single knots and flows 
(see e. g. 1994b). In order to calculate the inclination of each 
structure and to reduce the uncertainties, we have inferred from 
our optical data an ave rage v alue of both the tangential and ra- 
dial velocities. In Table lE.llFl the inclination and spatial velocity 
derived for each knot are reported (averaging both Ha and [S 11] 
values). A detailed analysis is given in Appendix[E] Inferred spa- 
tial velocities range from 25 to 125 km s"'. We obtained three 
different averaged values of inclination for the flows (58+3°, 
84+2°, and 67+3° for HH52, 53, and 54 flows, respectively), 
which support the proposed scenario of three distinct outflows. 
The outcome is illustrated in Fig. [12] where the measured knots 
are spatially plotted with their P.M.s, P.A.s, and inclinations (dif- 
ferent colours represent different range of inclinations with re- 
spect to the plane of the sky, a few more knots of HH 54 have 
been included combining [Fe 11] and optical data - see next para- 
graph). 

We also attempted to estimate the inclination angle for the 
structures encompassed by the [Fe 11] ISAAC slits, in order to 
define their proper outflow. To this aim we tentatively combined 
the [Fe 11] radial velocities with the atomic tangential velocities 
derived from the optical lines, keeping in mind that these species 
can trace different velocities and regions of the shock. The re- 
sults are reported in Fig. [12] as well. For knot Z in the HH54 
streamer we obtain an inclination angle of 72°+2°, in good 
agreement with the optical results. Both K and E (/=73°+2°) 
are also part of this flow. Knot Al could be part of the HH52 
flow (/=62°+4°), while the association of knot A3 is not clear 
(/=68°+8°). Also J, H2, and H3 are likely parts of this flow with 



Tab. IE. H is presented in Appendix |E] as online material only. 
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the inclination estimates of 58°+6°, 54° +5°, and 53°+3°, respec- 
tively. 

4. 7. Diagnostic results from spectroscopy 

From the ratio of the detected lines of our optical spectra, it is 
also possible to derive the physical parameters of the encom- 
passed knots. The selected optical transitions from [Oi], [Nii], 
and [S ii] were used to determine the electron density (A^e), elec- 
tron temperature (Tg), and ionisation fraction (Xg) of the shocked 
gas by adopting the so-called BE technique (e. g. BE99, Nisini 
et al. 120051 Podio et al. I2006I I. which in turn allows us to de- 
rive the total hydrogen density A^h (Ne/Xg) as well. The assumed 
chemical abundances are those from Esteban et al. (2004), ex- 
tensively discussed and also used in Podio et al. ( I2006I I. The re- 
sults are presented in Table[3](Col. 2-5). The observed lines have 
been dereddened using the optical extinction law of Cardelli et 
al. ( 119871 1 assuming an average Ay value of 2+1 mag for all the 
objects, as retrieved in previous works on HH54 from the NIR 
an alysis (Gredel |1994l Caratti o Garatti et al. |2006t Giannini et 
al. I2OO6I 1. The visual extinction in HH52 and 53 should have 
similar values as the extinction map in Porras et al. (120071 ) sug- 
gests (see their Figure 9). The errors of the physical parameters 
in Table [3] were derived by combining the uncertainties on both 
the line fluxes and extinction. 

The values appear quite homogeneous along and among the 
three flows, indicating the presence of very light and ionised 
jets, as sometimes observed in other HH flows, where shocks 
propagate in a medium of relatively low pre-shock density (see 
e. g. Podio et al. 120061 ). Actually, both ionisation fraction and 
electron temperature have high values, ranging between 0.2-0.6 
and 14 000-26 000 K, respectively. On the other hand, both hy- 
drogen and electron densities are extremely low with values of 
2x10^-2x10^ cm"-' and 50-790cm"^, respectively. Remarkably, 
such a low jet density for HH54 was predicted by the shock 
models of Giannini et al. (2006), supported by H2 and FIR 
observations. There, the emission in HH54 was interpreted 
as originated by a J-type shock with a magnetic precursor of 
A^H ~10^ cm"-*, which interacts with a relatively denser ambient 
medium (A^h ~10'* cm"^). 

4.8. Mass flux, mass, and luminosity 

The physical and kinematical parameters derived in the previous 
sections allow us to evaluate the mass flux rate of the knots (M) 
in the studied HH objects. In particular, together with the direc- 
tion retrieved from the P.M. analysis, M can help us to identify 
which exciting source is driving the outflow, comparing the mass 
flux rates of the outflow with the theoretical ejection rate of the 
candidate sources (also taking their evolutionary stage and bolo- 
metric luminosity into account). 

We used two different methods to infer M (for a detailed dis- 
cussion see Nisini et al. 120051 Podio et al. 120061 ). In the first 
method (A), the mass flux is estimated as: = fimnNHTrr^Vk, 
where jj is the average atomic weight, mu is the proton mass, 
A^H the hydrogen density, and Vk and rk are the velocity and 
radius of the knot, respectively. For this we have taken half of 
the length of the knot, derived from a 3 cr contour on the im- 
ages. Method A assumes that the knot is uniformly filled at the 
derived density, providing us with an upper limit of the mass 
flux. The second method (B) makes use of the luminosities of 
the optical lines ([Sii], [Oi], and [Nii] doublets), the hydrogen 
density, and the tangential velocity: - fj.mH{NHV)Vf/l,, with 



NhV - Lx(hvAifiY^) where Lx is the luminosity of the ele- 
ment X, for the selected transition, Ai and /; are the radiative rate 
and the fractional population of the upper level of the transition, 
Y is the ionisation fraction of the considered species with a total 
abundance of p| with respect to the hydrogen. Line fluxes were 
dereddened according to Sect. 14.71 The results from each doublet 
in each knot were averaged, obtaining a final value. The values 
of the mass flux derived from this method can be considered as 
lower limits (e. g. Nisini et al. 120051 Podio et al. 2006), because 
method B only measures gas that is sufficiently heated to radiate 
in the observed lines. A correction factor was applied when the 
slit width was less than the encompassed knot. The [S 11] images 
allowed us to derive such a quantity, however, for the other two 
atomic species. The dimension of the emitting region could not 
be inferred, and the same correction factor has been adopted, as 
well. In Table [3] (Cols. 7-8) the values of the mass flux derived 
from both methods are reported. The resulting mass fluxes are 
presented in Fig. \13\ as well, where for each flow, M of each 
group of knots is plotted as a function of the distance from the 
first knot of the flow. As expected M{A) is, on average, almost 
an order of magnitude greater than M(B) (in 2 out of 3 flows). 
Along the HH52 flow (bottom panel of Fig. [13]), M(A) osciUates 
between -2x10"^ and ~lO"^M0yr"' (between -3x10"^ and 
~3xlO"^M0 yr ' with method B). Mass flux along HH53 flow 
(M0(A) <lO"^M0yr"', Mq{B) <3x10"'' M© yr"' ) is more 
than an order of magnitude lower than in HH52. This could in- 
dicate a less massive and/or an older exciting source. The value 
in HH54X (along the streamer) is almost the same with both 
methods (-10"^ Mq yr"'). 

Using method A and considering an average hydrogen den- 
sity of lO-' cm"^ (see also Table [3j, we also estimated the mass 
flux for most of the other knot groups. The values are also re- 
ported in Fig. [13] as filled circles. M{A) has almost the same 
value in HH 52 E, F, and G (-2x10"^ Mq yr"') Along the HH 54 
streamer, the estimate of M{A) ranges between 3x10"' and 
~lO"^M0yr"' in Y and Z, respectively. On the main body of 
HH 54, we can evaluate a mass flux rate of about 9x, 4.5x, 3x, 
1x10"** Mq yr"' for A, E-K, H, and I, respectively. 

It is worth noting that such values are in excellent agree- 
ment with the mass flux rates estimated in Giannini et al. ( 2006] ) 
and with the linear momentum (P) obtained from CO in 
Knee (119921 ). considering the velocities and the mass flux rates 
measured in this paper This agreement between atomic and 
molecular momentum flux can only be explained if the jet is 
strongly under-dense with respect to the ambient medium (10- 
100 times). Then the efficiency of the momentum transfer from 
the jet (HH) to the medium (CO) is close to one. 

Moreover, to measure the energy budget of the flows, we 
computed mass and luminosity of the knots (for the detected 
species, i. e. HI, [S 11], and H2). To derive the mass, we assumed 
Mil = jjinnNnV, where fi is the average atomic weight, mu is the 
proton mass, A^h the hydrogen density, and V the volume of the 
knot. For those knots where the hydrogen density was not mea- 
sured, we assumed an A^h average value of lO^cm"^ (see also 
Table [3]). The luminosity of each element is L = Ancfil, where 
d is the distance and / the total intensity. The total intensity is 
derived from the dereddened line fluxes measured from narrow- 
band imaging. As in Sect. 14. 71 the adopted extinction value is 
Av=2+1 mag. For the [Sii] emission an NLTE model has been 
employed to derive the line ratio with respect to the sulphur dou- 
blet (see e. g Nisini et al. 2005 ; Caratti o Garatti et al. 2006 ). The 
input parameters are T^, Ne, X^, and Ay. Where Tg, N^, X^ have 
not been measured, average values of 20 000 K, 10"^ cm"^, and 
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Table 3. Measured physical parameters and mass fluxes of knots in HH 52, HH 53, and HH 54. 



knot 


Log (A',) 


X, (NJNn) Te 


Log (A'h) 


rk„„, M {Ay 


M (B)" 




(cm-3) 


K 


(cm-3) 


(10-'* Mgyr-') 


(10-'* Moyr-') 



HH52A1 

HH52A3-4 

HH52 D4 

HH52B 

HH52D2-D3 

HH53 A 

HH53B 

HH53C 

HH53C1 

HH54 C2 

HH54 C3 

HH54 G-GO 

HH54G1-G3 

HH54X3 

HH54X4A-X4B 



2.82±0.02 

2.2+0.2 

2.88+0.03 
2 45+0.18 

2.88+0.03 

2 60+""* 

2.46+0.11 

2 59+"- '3 

2.76+°;°? 

2.66+«- 

2.65+0.01 

2 00+" -' 

1.70+0.7 

2.6+0.3 



0.43+0.02 
0.35+0.06 
0.5+0.05 
0.40+0.03 

23+°-"'' 

"•^•'-0.06 

0.43+0.03 

31 +" "5 

"•-'^-0.09 

0.34+0.04 
0.32+0.09 

49+003 
0.50+0.07 
0.6+of 
35+0 0'' 

"•-'-^-0.11 

0.34+0.07 
0.45+0.05 



17500H 
23800H 
21000H 
21500H 
14200H 
21600H 
17300H 
18600H 
24200H 
22600H 
26000H 
19700H 
23000H 
25400H 
26000H 



:700 
:2600 

:1000 

900 

Mo 

■1800 
:800 
2500 
■1900 
1200 
■1100 
:3700 
2800 
■1600 
:7000 
800 
■1700 
:3000 
3600 
■3800 
:2000 



3.18+0.16 
2.9+0.5 
2.5+0.1 
3.28+°^^'^ 
3.094 



-0.76 

3.15+0.20 
2.92+0.03 
2.9+0.3 

3.06+0.20 

2.9+0.4 

2.85+°- 

2 17+o^5i 

^•^ '-0.48 

3.0+0.2 



1.5 
1.3 
1.6 
3.3 
2.4 
0.8 
0.8 
0.6 
0.6 
0.7 
0.7 
0.7 
0.7 
1.3 
2.1 



8.2+0.7 

1.5+0.4 

2.1+0.2 

16+2 

7+2 

1.1+0.2 

0.7+0.1 

0.3+0.03 

1.3+0.1 

0.3+0.05 

0.2+0.05 

0.6+0.03 

19+6 

0.3+0.08 

0.8+0.1 



1.1+0.1 

0.3+0.03 

0.2+0.03 

3+0.3 

1+0.1 

0.2+0.02 

0.3+0.02 

0.07+0.01 

0.3+0.06 

0.3+0.05 

0.1+0.05 

0.05+0.01 

3+0.5 

0.4+0.2 

0.6+0.2 



Notes: 

" M derived with method A (see text). 

* M measured from [Sn], [Oi], and [Nii] line luminosities (average value) multiplied by a correction factor that takes the entire area of the knot 
into account (method B, see text). 



4.9. The candidate driving sources 
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Fig. 13. Mass fluxes derived for HH52 (bottom panel), HH53 
(central panel), and HH 54 flows (upper panel). Circles and tri- 
angles refer to M values inferred from method A and B, respec- 
tively. Filled circles refer to M estimates obtained with method 
A, assuming an Nh average value of lO-' cm"-'. 



0.4, have been adopted, respectively. The HI intensity was ob- 
tained from the Ha line, assuming that the radiated energy is 
emitted under Case B recombination, and using emissivity val- 
ues from Storey & Hummer ( 1995 ). Finally, for the molecular 
hydrogen, when observed, we consider an LTE gas at T-2 000 K 
(see e. g Caratti o Garatti et al. |2006l l. 



As seen in Sect. 14. II our H2 maps do not reveal any molecular jet 
close to the possible driving sources, therefore we do not have 
any direct indication of which are the driving sources of these 
flows. We thus selected the possible candidates on the basis of 
several YSO catalogues (i. e. IRAS, 2MASS, and MSX) and pub- 
lications (see e. g. Young et al. |2005i: Spezzi et al. 2007; Alcala et 
al.[2008 ; Spezzi et al. 2008 ), on the IR colours of the YSO in the 
Spitzer images, and on their positions with respect to the HH ob- 
jects. In Fig.|2l the 24fj.m Spitzer-MIPS map of the region indi- 
cates the selected candidates and the location of the studied HH 
objects. There are five IRAS sources in the region with a position 
barely compatible with the derived P.A. of the HH flows. These 
young sources have been deeply analysed in the literature (see 
i. e. Spezzi et al. |2007l Alcala et al. 12008', Spezzi et al. |20U8] l. 
where the authors report their main characteristics and spectral 
energy distributions (SEDs). From the Spitzer IR colours, we 
recognised two more possible candidates that have been added 
to our list. We note, however, that these two objects were not in- 
cluded in previous Spitzer studies of the region, because they do 
not fit some of the selection criteria used (see Spezzi et al. 2007|. 
In particular the SED analysis in the C2D catalogue indicates 
that Source 1 (see Fig. |2]l could be a galaxy. The sources, to- 
gether with short comments, are reported below. The SEDs of 
the eligible candidates are shown in Fig. [141 They have been ob- 
tained combining 2MASS, Spitzer, IRAS, and literature data. 

IRAS 12496 - 7650 (DK Cha) is a Herbig Ae star 
with a high velocity and weak CO outflow (RA.~10°-15°, 
Knee. 1 1992"]) . appar ently not related with HH52, 53, or 54 out- 
flows (Knee, !1992| l. even if its estimated mass accretion rate 
{Mace ~10-*'Moyr-\ Alcala et al. 2008) woul d mak e it an ex- 
cellent candidate (see discussion in Alcala et al. 120081) . Our RM. 
measurements of the three flows indicate, however, that this can- 
not be the driving source of any of these flows. 

IRAS 12448 - 7650 is a Clas s III YSO with a Mace ~10'^- 
10-'' Mq yr-' (Alcala et al. 12008] ) and located -68° with respect 
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to HH52 flow, i. e. 10°-15° of misalignment. Moreover, its opti- 
cal spectrum (Alcala and Spezzi private communication) shows 
no Li I absorption feature, indicating that the source has an age 
>2xl0'' yr (see e. g. Magazzu et al |1992| l. 

IRAS 1 2416 - 7703, a Class II YSO (Young et al. l2005l 
Alcala et al. l2008l l. has the correct position angle (~56°) to drive 
the HH52 flow. The SED analysis reveals that the source has 
a bolometric luminosity of about 8 L0 (see Fig. [T4l i (a simi- 
lar value is by Spezzi et al. I2008I I. In theory, this would be 
the most reliable candidate for the HH52 flow. Nonetheless, 
its optical spectrum (Alcala and Spezzi private communication) 
has a weak Ha emission but no Lii absorption feature, thus 
IRAS 12416 - 7703 seems too old to be the driving source. It 
is puzzling, however, that no other YSOs are present in this re- 
gion (see Figs.|2]and 9 in Alcala et al. l2008l l. 

IRAS 12500-7658 is a Class I (Che n et al.[T997| source with 
an Lboi of 0.5 L0 (Young et al. 2005) with an estimated Mace 
of 2x10 ''Moyr^' (Alcala et al. ,2008i) . Indeed this is the best 
candidate for the HH54 flow (M ~lO~^M0yr"', see Fig. [13) 
because of the P.A. of ~20° with respect to the HH object and 
the MjMacc (0.1-0.01) ratio, compatible with a Class I source. 

IRAS F12488 - 7658/C13, located two ai'cmin utes w est of 
the previous source, is a Class III YSO ( Voung et al. 1200 11 1 with 
a M„cf4xl0 yr ' (Alcala et al. l2008] l and has a P.A.~30° with 
respect to HH 54. The misalignment, the age, and the Mace of the 
source seem to indicate that this is not the driving source. 

Finally, we detect in the Spitzer images an MIR source al- 
most coincident with HH53 (indicated as Source 2 in Fig. |2] 
a = 12''55'"17* and 6 =-76°57'29'.'3), appearing as a faint stai- in 
our optical, H2, and 2MASS images. A faint emission is also vis- 
ible in the IRAC images (partially detached with respect to the 
HH53 knots), while this object is not resolved in the MIPS im- 
ages and part or all the measured flux could originate in the HH 
itself. The coordinates of this object are almost the same as the 
infrared source HH53*1 reported by Sandellet al. ( 1987) and in- 
dicated as the possible driving source of HH 53. On the contrary, 
Graham & Hartigan d 19881 1 firmly rejected this hypothesis be- 
cause their observations showed a blueshifted emission both to 
the west and to the east of the star (we know, however, that such 
emission could originate in the HH 52 flow) and because it did 
not show any unusual strong emission feature in the star spec- 
trum. Furthermore, the optical photometry of Spezzi and Alcala 
(private communication) indicates that the object would not be- 
long to Cha II and could be a background star. Nonetheless, our 
SED analysis (Fig. [141 from the Spitzer photometry and un- 
der the uncertain hypothesis that the MIR comes entirely from 
the object) would indicate that it is a Class II YSO, with a 
Lboi ~10"^ Lq. We tentatively speculate that this could be the ex- 
citing source of HH 53. The mass ejection rate measured on the 
knots (3x10"^-10"^Mq yr"') appears, however, too high com- 
pared to the Class and the luminosity of the presumed source. 
Indeed a more in-depth spectroscopical analysis is needed to 
solve this conundrum. Nevertheless, it is worth noting that no 
other candidates are detected east of HH53. In conclusion, two 
of the three identified flows (i. e. HH52 and 53) still have not a 
reliable driving source candidate. 

5. Discussion 

5.1. Variability 

5.1 .1 . Flux and velocity variability of knots 

Our analysis indicates that flux and velocity variability is ob- 
served in knots along the flows mostly in coincidence with in- 
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Fig. 14. SEDs of the discussed candidate driving sources. 



teracting knots or working surfaces, where pairs of knots inter- 
act by accelerating and decelerating the gas in front and behind, 
respectively, or where the jet/streamer hits the slow moving ma- 
terial ahead. On the other hand, knots with a constant or slightly 
variable flux do not show relevant changes in velocity. It is there- 
fore reasonable to investigate if and how such a flux variability 
is connected to the observed velocity variations. We can argue 
that part of the kinetic energy of the impacting knot is trans- 
ferred to the preceding knot or to the surrounding medium, and 
then converted in radiated energy. However, it is worth noting 
that here we do not consider how the energy is transformed and 
which processes are involved. In a forthcoming paper (De CoUe 
& Caratti o Garatti, in prep.), we report how the observed vari- 
ability can be reproduced by MHD simulations, when knots in- 
teract along the flow. 

In order to check how much of the kinetic energy is radiated 
during this process, we therefore computed the energy budget for 
those knots that are interacting and that clearly show both flux 
increment and velocity decrement. We estimate both the amount 
of kinetic energy lost in the interaction (AK) and the variation in 
the radiated energy during the time interval of our observations 
(AErad)- The kinetic energy lost by the knot can be expressed as 



AKk = 1 llMkAv' 



(1) 



where Mk is the mass of the knot and Av^ is the diflFerence be- 
tween the initial and final squared velocities. 
On the other hand, AErad can be written as, 



AE,. 



'■ad 



ALdt) 



(2) 
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, which is the integral of the luminosity variations over the time 
interval of observations, summed for all the detected species (i. e. 
HI, [S ii], and H2). Knot luminosities, and masses were estimated 
as described in Sect. 14.81 The results are reported in Table |4] 
where h.Erad of each element, and the total t^Erad are indicated 
for each analysed knot A/T. 

As a result A^rarf is always smaller than A/T. This is pos- 
sibly due to the longer cooling time of the gas with respect 
to the interval of our observations. Moreover this is also be- 
cause we are not taking into account the l^Erad contribution from 
other species like Fell, 01, or CO, for example. We can give a 
rough upper limit estimate of the Fell contribution, if any, for 
example. The brightest [Fen] emission in HH 52-54 complex 
is observed in HH54B, where a flux of ~8xlO '^ergs ' cm"^ 
(1.64 yum line) has been measured (see Gredel 1994 Giannini et 
al.[2006 Caratti o Garatti l2(306l l. This corresponds to a luminos- 
ity of Lfeii -4x10^^ erg s"' (adopting the same physical param- 
eters used for the other ions). This means that even assuming a 
luminosity variability equal to the knot luminosity the radiated 
energy in Fell would be 10^^-10^^ erg, which is one or two or- 
ders of magnitude less than the HI contribution. The contribution 
from [S 11] is usually smaller, down to two orders of magnitude 
(see TablelDl. Therefore, among the considered atomic lines, HI 
radiates most of the energy observed in such a variability. This 
is well visible from the flux variability in Figures [8]|9] (see also 
Figures in Appendix ICTi. where Hot usually exhibits the largest 
flux variability. This is possibly due to the different cooling times 
of the two atomic species (a few years for HI and ~30 years for 
[Sii], see e.g. Hartigan et al. 120011 12005I I. As a consequence, 
usually Ha traces morphological changes better and faster than 
[S 11] does. 

Conversely, the H2 emission, when observed, partially con- 
tributes to the gas cooling, but with longer timescales. We note, 
however, that only part of the H2 variability was estimated (i .e 
form the warm component ~2000 K traced by the 2.12 fim line), 
and we are likely missing the contribution of the cold H2 compo- 
nent that can contribute up to 50% to the whole radiative cooling 
(Caratti o Garatti et al. 120081 ). 

It is not possible to properly evaluate the remaining molecu- 
lar contribution, as for the CO. However it should not be greater 
than the atomic one, since the observed medium is partially 
ionised. Therefore it is reasonable to conclude that a large part 
of the kinetic energy lost in between two working surfaces is 
converted in radiated energy. This also means that the energy is 
rapidly converted, the shock dissociates the H2, and we observe 
the HI emission. After a few years, the temperature drops, and 
the H2 emission is then detected. 

However, in some cases, the kinetic energy is not completely 
radiated but can be partially transferred to the medium or to the 
knots ahead of the flow, as shown in the interaction between 
knots HH54G0 and G, for example. Here, the kinetic energy 
lost by the bullet (GO) is mostly ceded to the target (G) and only 
partially lost (~20%) in the radiative process. In fact, the energy 
balance of the system {hK{GQ)-AK{G)-hEradiGQ + G)) is close 
to zero, as can be derived from Table |4] 

Finally, there are some knots that show flux variability but 
no change in velocity (inside the error bars). In some cases, 
where the flux decreases, we are observing the radiative cool- 
ing of the knot, but the momentum of the knot is conserved. In 
other cases, where the flux increases (see e. g. HH 54 A 1 and A8 
in AppendixICli, we have estimated the amount of kinetic energy 
needed to produce the observed l^Erad- As a result the variation 
in velocity that produces such a AK is well below the error bars 
that we measure in our P.M.s. 



5.1 .2. Cooling time and velocity variability 

Our set of observations indicate apparent acceleration and de- 
celeration for several knots, but just a few of them seem to be 
genuine. Often the observed velocity changes do not reflect true 
fluid velocity variations, but merely phase velocity effects (see 
e. g. Bally et al. I2002I I. wherein some portions of the post-shock 
gas fade and others brighten. This is particularly true when the 
elapsed time between two examined images is longer than the 
post-shock cooling time of the analysed knot. 

We can thus compare the timescale of the observed variabil- 
ity with the expected atomic and molecular cooling times. The 
cooling time is 



tc. = 



nk^T 

\dT/dt\ ~ A 



(3) 



where T is the temperature, the Boltzmann constant, and A 
the cooling rate (in erg cm"^ s"' H To estimate t^, we use the H2 
cooling rates from Shull & HoUenbach (1978), and the atomic 
cooling coefficients from Cox & Tucker (1969) and Gnat & 
Sternberg (120071) . For the H2 we consider a range of densities 
(from 10^-10^ cm-3) and temperatures (from 2000 to 3000 K). 
As a result, the H2 cooling time ranges from ~1 to 9yr, where 
the shorter value is obtained for higher temperatures and densi- 
ties. For the atomic component, we consider an average A^h of 
lO^cm-3 and ranging from 10000 to 25 000 K. We thus ob- 
tain 1< tciHa) >6yr. The cooling time of the [S 11] is almost one 
order of magnitude longer than t^iHa). 

The complex proper motion variations may then indicate the 
combined effects of photometric and morphological changes, oc- 
curring on timescales comparable to the cooling time in the post- 
shock gas, combined with the detection of collisions of fast, frag- 
mented ejecta overrunning slower moving fragmented ejecta. 
The highly fragmented morphology of these shocks and the very 
short acceleration/deceleration timescales (-lOyrs) indicate that 
various types of instabilities (cooling instabilities, Vinsniac in- 
stabilities, etc) may have caused the flow elements to fragment 
into tiny, dense clumps surrounded by lower density plasma or 
gas. As a fast shrapnel over-takes slower shrapnel, both acceler- 
ating and decelerating post-shock fluid elements are observed. 
Such timescales could be actually good measurements of the 
size-scales of the clumps in the colliding fluids. 

5.1 .3. Direction variability and flow deflection 

Our P.M. analysis clearly indicates that part of the HH 52 flow is 
deflected. Some knots in the HH 54 G group have been observed 
to change their direction 'in real time', due to the interaction 
with other colliding knots moving faster along the flow. In this 
case, most of the kinetic energy of the impacting knot (bullet) is 
transferred to the knot ahead (target), as shown in the previous 
section. To cause the deflection, the bullet trajectory is slightly 
off-axis with respect to the target. All this implies a time variabil- 
ity in the ejection velocity and a small precession of the flow, as 
weU (see e. g. Raga et al. [T993l Volker et al |T999] l. 

Other knots, such as HH 54 A, C, M, or HH 52 G, F exhibit- 
ing different P.A.s of 20°-40° with respect to the main flow, 
were probably deflected earlier, and it is not evident how this 
bending originated. However, it is not clear whether this 'bullet- 
target' mechanism would remarkably deflect massive features 
like HH54A. Such a deflection is more likely the result of a 



' ' alternatively the cooling coefficient A,, (in erg cm^^ s ' ) can be used 
and Eq.[3]becomes: tc ~ ^ 
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knot hK AE,„d{Hl) AE,-„d(SH) AE,ad(Hi) AE,-aAtot) 



ID 


(103>*erg) 


(10^^'* erg) 


(10'** erg) 


(lO'** erg) 


(10^** erg 


HH52B1 


70 


1.1 


0.08 


0.2-0.4 


1.3-1.5 


HH52C1 


17 


0.8 


0.18 


04 


1.3 


HH52C2 


3.5 


0.35 


0.1 




0.4 


HH54H1 


8.8 


2.9 


0.2 




3.0 


HH54H2 


9.6 


2.8 


0.12 


0.9 


3.8 


HH54H3 


57 


2.7 


0.2 


2 


4.8 


HH54G 


-18" 


4.0 


0.02 




4.0 


HH54G0 


26 


0.8 


<0.01 




0.8 



Notes: 

" here the kinetic energy is increased from the interaction with GO. 



collision of the flow with a dense cloud or dense clumps in 
the medium, as also observed in HH 110 (Reipurth et al l 19961 1 
and in several numerical simulations (see e. g. de Gouveia Dal 
Pino [19991 Raga et al. [20021 Back et al. I^OSb . In this sce- 
nario the deflected beam initially describes a C-shaped trajectory 
around the curved jet/cloud contact discontinuity; later, when the 
jet has penetrated most of the cloud/clump extension, the flow 
resumes its original direction of propagation. This could be the 
reason we observe both a deflected (HH 54 A) beam and a con- 
tinuation of the flow (HH 54 G and C). 

The observed velocities in the deflected beam are between 
30% and 60% lower than those detected in HH52. This is con- 
sistent with the values predicted by the simulations (see e. g. de 
Gouveia Dal Pino [19^91 Back et al. l2008b with v/,„=v,„cos0, 
where 9 is the deflection angle, v/,„ and v,„ are the deflected 
and incident velocities, respectively. Also the predicted ratio be- 
tween jet and ambient mediu m (»/ /»m <10"^) agrees with our 
observations (see Sects. [4!7l and l4.8l l. Therefore the chaotic struc- 
ture and PM.s measured in HH 54 should be interpreted not only 
as the simple superposition of two distinct outflows but also as 
the result of a deflection in the HH 52 flow. 



5.2. Length, dynamical age of the outflows, counterflows 

Assuming that //?A5 12416-7703, Source 2, and 12500- 
7658 are the driving sources of HH52, HH53, and HH54, re- 
spectively, we can infer two important parameters of the flows 
from our kinematical analysis: the spatial length and the dynam- 
ical age (Tdyn), which is the ratio between the distance of the 
knot(s) from the YSO (corrected for the inclination angle /) and 
the total velocity. The inferred results are reported in Table |5] 
Additionally, Fig.[T5]summarises the results, showing the candi- 
date exciting sources, the outflows, and their orientations. Two 
out of the three outflows, namely HH 52 and HH 54, would be 
parsec scale outflows, with a length of about 4 pc and 2.2 pc, re- 
spectively. The dynamical age of the flows would then be up to 
6x10^* yr (HH 52 flow, knots HH 54 A and C), and up to 4x10* yr 
(HH54 flow, knots E, K). On the other hand, if we assume 
Source 2 as the HH53 driving source (but this is highly ques- 
tionable), here the flow has a smaller extension 0.1 -0.2 pc and a 
shorter dynamical age (<10-' yr). 

The detection of parsec scale outflows from Class & I 
(e. g. Reipurth et al.[ IMl Eisloffel [2000b as from Class II YSOs 
(McGroarty & Ray I2004l l is quite common. It is also not sur- 
prising that these flows have 4-6x10"^ years of age. What is in- 



Table 5. Lengths and dynamical age of the flows. 



knot 


Possible 


Length 


'^dyn 


ID 


Exciting Source 


(pc) 


(lO'yrs) 


HH52A 


IRAS 12416-7703 


3.9 


39 


HH 54 C 


IRAS 12416-7703 


4.2 


63 


HH54G 


IRAS 12416-7703 


4.2 


41 


HH54A 


IRAS 12416-7703 


4.2 


63 


HH53 A-C 


Source 2 


0.1-0.2 


0.3-0.8 


HH54 Y 


IRAS 12500-7658 


2.1 


18 


HH54 K-E 


IRAS 12500-7658 


2.2 


43 



deed peculiar is that no further emission is detected between the 
parsec scale flows and the possible exciting sources. This could 
mean that the ejection (and therefore also the accretion) activ- 
ity in two of the YSOs suddenly stopped or at least was reduced 
several thousand years ago. A different explanation could be that 
part of the flow is not emitting (i. e. it is neutral), and it only be- 
comes visible when it shocks the slow material moving ahead. 
This seems to be the case of the HH 52 B and C regions, to the 
rear of the HH 52 flow, where we observe the formation of new 
knots. 

Another puzzling question is the absence of a counterjet in 
all three flows. The red-shifted CO emission between HH 52 and 
HH54 (see Knee [T992l l could be the HH53 missing counter- 
flow. However, we do not detect any jet in our images. Also for 
the two parsec scale flows, no emission is observed along the 
possible counterflow axes in Spitzer, 2MASS, or DSS images. 
Alcala et al. (2008) report the discovery of a faint Ha emis- 
sion, located -20' to the SW of IRAS 12496 - 7650, but it is not 
aligned to the HH 54 flow. Finally, the Spitzer maps do not cover 
the hypothetical coordinates of HH52 and HH54 counterflows. 
According to CO maps of the Cha II cloud (see e. g. Hayakawa et 
al. 2001 ), both these coordinates are close to the southern cloud 
edge. Thus, because of the very low medium density, it is possi- 
ble that the ejected matter cannot produce any shock. 

6. Conclusions 

We have presented a detailed multi-epoch (20 years baseline) 
kinematical investigation of HH52, 53, and 54 at optical and 
near-IR wavelengths, along with medium (optical) and high- 
resolution (NIR) spectroscopical analysis, probing the kinemat- 
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u 

Id 



20;«f0.0 




THAS 135C0 765* 



MIPS 24 



5<:«,0 54;O0,O 52iW.O I2;»:0a-9 J8;00-0 4*;00.« 4l;(»J 

RJL. (J2000) 

Fig. 15. Summary figure. Spitzer-MIPS 24 jjm map of the investigated Cha II region, indicating the candidate exciting sources, the 
outflows, and their orientations. A rectangular box indicates the position of the upper-right inset, where the Ha EMMI 2006 image 
along with the H2 SofI 199 image contours of the three outflows are reported. 



ical and physical conditions of the gas along the flows. We in- 
vestigated flux and velocity variability with time of the flows, 
reporting the first detections of acceleration and deceleration of 
protostellar jets. We find that there is a correlation, at least in 
some cases, between the flux and the velocity variabilities ob- 
served along the flows. Finally, we discussed the possible flow 
exciting sources. The main results of this work are the following: 

- The P.M. analysis reveals that there are three distinct out- 
flows, partially overlapping. The first outflow (HH54 flow) 
is oriented NNE with a position angle of ~22°, delineated by 
the HH 54 streamer, and includes some of the knots in the 
main body of HH 54. The second (HH 52 flow) follows an 
NE direction with a position angle of ~55°, grouping HH 52 
and part of HH 53 and 54 knots. The third (HH 53 flow), de- 
lineated by the three brightest knots of HH53, is moving in 
an E-W direction. The HH 52 and 53 outflows partially over- 
lap and seem to converge on the HH 54 bow shock. 

- With a few exceptions, the derived P.M.s have similar val- 
ues in the the optical filters (Ha, [S 11]) ranging between 
0.007 "yr"' and ~0.097 "yr"', corresponding to a tangential 
velocity between 6 and 82kms"' at a distance of 180pc. 
On average the H2 P.M.s have lower values (from 0.007 to 
0.041,"yr"'). 

- Measured radial velocities of the atomic components (both 
optical and NIR) range between -20 and -120 kms ' with 
variations similar to the P.M.s. Spatial velocities of the knots 
vary from 50 kms"' to 120 kms"'. The inclinations of the 
three flows are 58+3°, 84+2°, and 67+3° (HH52, HH53, 
and HH 54, respectively). 

- In 20 years, about 60% of the observed knots in the flows 
show flux variability. Moreover, our set of observations indi- 



cates acceleration and deceleration in several measurements, 
but just some of them seem to be reliable. These knots are 
working surfaces or interacting knots. In this case a rele- 
vant flux variability is observed as well. We argue that both 
flux and velocity variability are related and that part of or all 
the kinetic energy lost by the decelerating knots is succes- 
sively radiated. However, the complex proper motion varia- 
tions may also indicate the combined effects of photometric 
and morphological changes, occurring on timescales compa- 
rable to the cooling time in the post-shock gas, combined 
with the detection of collisions of fast, fragmented ejecta 
overrunning slower moving fragmented ejecta. 

- The physical parameters from the diagnostic are quite ho- 
mogeneous along and among the three flows. The analysis 
indicates the presence of very light (A^h ~10^ cm"^), ionised 
(Xe -0.2-0.6), and hot (T^ ~ 14 000-26 000 K) jets, impact- 
ing against a denser medium. 

- The measured mass flux rates of the knots span from 
10"^ Mq yr"' to some 10"^ Mq yr"', where the highest val- 
ues are detected in the HH 52 flow. 

- Several knots are deflected, at least in the HH52 flow. At 
least for a couple of them (HH 54 G and GO), the deflection 
originates in the collision of the two (what we call a bullet- 
target deflection mechanism). However, for the more mas- 
sive parts of the flow, the deflection is more likely the result 
of the flow collision with a denser cloud or clump. 

- We investigated the possible driving sources of the flows, 
comparing knot P.M.s and mass flux rates with the position. 
Class and bolometric luminosity of the source candidates. 
We indicated three possible candidates, IRAS 12416 - 7703, 
Source 2, and IRAS 12500 - 7658 for HH52, HH53, and 
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HH 54, respectively. Only IRAS 12500 - 7658 is, however, a 
reliable candidate. 
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Appendix A: Morphology 

A.1. HH52-53 

Our high resolution images allow us to disentangle the struc- 
ture of the studied HH objects. Figure H] (right panel) and lA.ll 
show the HH52 and 53 region in the three different filters. In 
Ha and [S ii] HH 52 exhibits a roughly bow shape, where knots 
labelled A can be identified as the head of the shock (nucleus). 
The brightest knot (Al) is surrounded by several smaller struc- 
tures, with knots A3 and A4 placed behind, while A5 and A6 
are in front, elongated towards a group of separated knots (Dl, 
D2, and D3) roughly E-W aligned. A faint emission (knot D4) 
is placed just in front of the nucleus. Groups of knots B and 
C constitute the right and left wings of the bow, respectively. 
The right wing appears particularly fragmented with knots Bl- 
B3 representing a first condensation close to the nucleus. B4-B9 
are located behind. 

The H2 emission (Fig. lA.ll bottom) shows a more amor- 
phous shape. A semi-ring shaped region, roughly coincident 
with knots A2, A3, and A4, is located just behind the nucleus. 
Three more bright spots identify the two wings. The first emis- 
sion coincides with CI on the left wing, the second, on the right 
wing, extends to an area corresponding to B2 and B3 in the opti- 
cal, and the third to B5 and B6 on the right wing as well. Finally, 
fainter and diffuse emission is detected on the Dl, D2 and D3 
knots. 

In addition to the three HH53 knots already known (A, B, 
and C), we detect several more emissions coming from this re- 
gion. Three more knots, mainly emitting in the Ha line (Fig.|4|, 
are found along the main flow, forming an elongated S-shape 
chain of E-W knots. We label CI the emission between C and 
A, and Bl and B2 those close to B. These objects are not de- 
tected in our H2 images and two of them (knots CI and Bl) are 
just barely visible on the [S 11] images. More diffuse emission, 
detected in all three filters, appears superimposed on the HH53 
main flow, connected to the HH52 streamer. HH53I lies SW, 
towards HH 52, while H and E are situated NE in the opposite 
direction. Two brighter structures are located farther NE along 
the streamer Group F, with a jet-like appearance in the optical 
filters, is approximately elongated towards NE, also visible in 
H2. Knot G, only visible in the optical, is approximatively lo- 
cated at the origin of the continuous emission of the streamer 

A.2. HH54 

In Figures 13] (left panel) and lA.2l we show the HH 54 images and 
contours in Ha, [S 11] (left panel of Fig. |A.2| l and H2 (right panel 
of Fig. IA.2b . The morphology of the object appears extremely 
complex. Our images reveal several substructures. For the sake 
of simplicity, we divide the description of the region in two parts, 
the streamer, composed of groups X, Y, and Z, and the main body 
of HH 54, made up of the remaining knots. 

HH 54 streamer shows an uncommon shape of a double he- 
lix, clearly visible in Ha. The brightest helix is also well delin- 
eated in [S 11], appearing as a wiggling jet fragmented in sev- 
eral knots. The morphology of the streamer closely resembles 
that of the HH 46/47 outflow (see e. g. Eisloffel & Mundt 1994b , 
Heathcote et al [T996) , except for the presence of a second faint 
helix. The base of the streamer is XO, with knots X2 and X3 de- 
lineating the first bend of the bright helix. The path follows with 
X4. Yl and Y2 delineate the second bend of the jet, that con- 
tinues with Y4. The faint helix is delineated by XI, Y3, Y8, Y7 
and Y6. Both the bright and faint structures end in two shocks 



(Z and Zl, respectively), that barely appear as a bow shock fea- 
ture in H2 (Fig. |A.2l right panel). Molecular hydrogen emission 
is also detected in knots XO, X2, Yl, and Y2. A further knot (T) 
complementary to the ionic emission is observed along the jet. 
Finally, some filaments (X4 A, Y5, R), detected at optical wave- 
lengths, are located aside, detached from the main flow. 

The morphology of the main body appears even more com- 
plicated than the streamer. In addition to the original knots (from 
A to K), observed in previous papers (Schwartz 1 19771 Schwartz 
& Dopita 1980, Sandefl et al.'1987), we detect more features (la- 
belled L to S). Moreover, several knots present sub-structures, 
labelled here with numbers. 

Group A at the bottom appears as bow-shaped in the optical, 
with a P.A. of about 75°. In H2 the group exhibits emission only 
in the northern part of the structure, only partially overlapping 
the atomic emission. A fainter arc-shaped emission (group Q), 
mainly visible in Ho-, is located east. Knots P (detected in the 
optical) and S (detected in H2) are located west, and together 
with knots O (farther north, visible in the optical) form the far- 
thest point of the HH 52 streamer. 

Towards north, more groups are detected. From east to west, 
almost coincident in all three wavelengths, we detect groups K 
and E, slightly elongated northward, then group B, and, finally 
group F, detected only in the optical. Farther NE an arc-like 
structure is identified by groups J and M in both optical and NIR. 

In the north-eastern region, knot I and group H appear as 
optical jets with RA.s of ~45° and ~20°, respectively. An arc 
emission in H2 is delineated ahead of I. 

In the north-western region, we detect group G, composed 
of five subsequent knots only detected in the optical (mainly in 
the Ho- filter, see also Fig.|4]i, with a P.A. of ~55°. At the end of 
this flow another group of knots (C) is observed, with CI and C2 
already known, and the new knot C3 located on a straight Une 
connecting knots G and C2. 

It is worth to note that the identification in the NIR of knots 
C accomplished by Gredel (1994) (and then used in Caratti o 
Garatti et al. 2006 and in Giannini et al. 2006) is erroneous. 
Indeed the knots that were labelled in H2 as CI and C2 have a 
similar appearance in the optical, but are situated about 6" SW. 
As a consequence, the knot previously known as C 1 in the NIR 
coincides with knot H2 in the optical, while C2 roughly corre- 
sponds to CI in the optical (see Fig.|4land |A.2t . In this paper we 
use this nomenclature to avoid further confusion. 



Appendix B: Proper motions 

B.1. HH52-53 

Knot position angles along HH 52 flow have a quite wide spread 
in their values (30°-80°). The bulk of HH52 and other knots, 
as HH 53 CI, I, and H, have a PA. of -50°. On the other hand 
HH52A3, A4 and a few more knots as HH53F and G, pro- 
ceed towards HH 54 with different P. A.s between 60° and 80°. A 
completely different behaviour is found for the features located 
on HH 52 wings in the rear of the flow. Although following the 
bulk of the flow, they show a high degree of variability during 
the 20 years of observations, as reported in Appendix ICl Such 
a variability, both due to local instabilities inside the flow and 
fast moving shocks overtaking HH 52 from behind, makes their 
identification in different epochs and a correct evaluation of the 
PM.s and PA.s sometimes difficult (see also Fig|5]and|6]l. 

In Fig. IB.ll H2 proper motions are reported (see also Tab. IB. 31 
and lB.4b . Significant measurements in HH 52 were possible only 



A.Caratti o Garatti et al.: Time-variability in the Clia II outflows, Online Material p 3 




Fig.A.l. HH52 and HH53 regions with [Sii] (EMMI 2006) (top), and H2 (Sofl 1999) (bottom) fihers. The labels indicate the 
position of the knots, including the newly detected ones. The contour levels of the [S 11] are 3, 10, 20, 30, 50, 60, 80, 90, 100, 150, 
200x cr (~ 4 X 10"'^ erg s^' cm^^ arcsec^^). H2 levels are 3, 10, 20, 30, 40, 50, 60x cr (-10"'^ erg s^' cm"^ ai-csec^^). 



for the three bright spots located in the wings. B2-B3 and CI fol- 
low the flow as in the optical, with similar tangential velocities 
and position angles. 

HH 53 A, B and C, in the HH 53 outflow, exhibit small P.M.s 
with P. A.s around 270° also in the NIR. Here, however, the mea- 
surements are less significant due to the larger errors. 

B.2. HH54 

Proper motions in the HH 54 streamer range between -0.01 and 
0.07" yr"' (i. e. Vtan 8 and 61 kms"'), where the highest veloci- 
ties (between ~30 and 60kms"') are detected in the inner part 
of the flow and exhibit alternating P.A.s (322°-69°), likely due to 
precession, with directions that oscillate around an average value 
of -20°. 

Some features, located on the edge of the streamer, are not 
collimated and are drifting apart (P.A.~90°) at a lower veloc- 



ity. At the tip of the streamer group Z appears almost stationary 
in Ha (vtan ~10kms"'), while in [Sii] has a higher value of 
~30kms-'. 

Moving towards NE we observe a decrease in the tangen- 
tial velocities from groups Q (~40kms"') to K and E (~20 
and 15kms"', respectively). From the P.A.s and inclinations 
(see Sect. 14.61 ) we can reasonably associate these knots with the 
HH54 flow. 

Group A appears as an expanding bow shock, towards ENE 
with an average direction of about 80° and a mean Vtan of 
~30kms"'. Knots Al, A2, A6 and A8 exhibit some variability, 
as reported in Appendix ICl Because of the morphology, the di- 
rection, and the average inclination (see Sect. 14. 6t . the structure 
appears to be connected to knot P and to the HH 52 flow. Indeed 
part of the flow along the HH 52 streamer curves and shocks this 
region, producing group O as well (see also Fig. [TJ. H2 emis- 
sion in HH54 A (Fig. lA.2t appears as an expanding arc, and its 
motion (vtan ~30kms"') has a different orientation depending 
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Fig. A.2. HH54 regions with [S ii] (EMMI 2006) (left) and H2 (Sofl 1999) (right) filters. The labels indicate the position of the 
knots, including the newly detected ones. [S n] contour levels are 3, 10, 20, 30, 40, 50, 60, 70, 80, 90, lOOx the standard deviation 
to the mean background (~ 4 x 10"'^ erg s^' cm~^ arcsec"^). H2 contour levels are 3, 10, 20, 30, 40, 50, 60x the standard deviation 
to the mean background (~10""' erg s"' cm"^ arcsec"^). 



on the position along the structure (Fig. IB. 11 1. Proceeding from 
east to west, we separately measured four regions, indicated in 
Tab. IB. 41 as Ai, A2, A3, and A4. The directions of the motions 
have different P.A.s with respect to the optical counterparts and 
it is not clear if the molecular emission comes from material 
swept up by the bow shock on one side (see e. g. HH 219, Caratti 
o Garatti et al. I2004I I or from the HH 54 streamer 

The faintest component of the HH52 streamer follows a 
straight trajectory (of ~55°) and impacts ahead, producing 
groups HH 54 G and then C. We observe a sudden deceleration 
of the flow along the HH52 streamer (from ~65 to 30kms ') 
as it collides with the leading material, that is partially deflected 
sideways (P.A.~80°). GO is part of such a fast flow, that shocks 
at the end of a slow moving component (knot G), changing its 
direction (see also Sect. 14.41 ). This dynamic is also visible in 
knots Gl and G3, that are not aligned with the flow, moving 
with a PA. of ~80° and ~65°, respectively. Group C ahead of 
G is proceeding in the same direction at similar tangential ve- 
locities (25-30kms"'). Surprisingly, knot C3, at the tip of the 
flow, slightly moves in the opposite direction in both optical and 
NIR filters (PA.~250° and Vtan ~15kms-'). This behaviour is 
observed also in knot Kl, on the opposite side of the HH object, 
showing nearly a reverse motion with respect to the main flows. 

Velocities of knots M are quite low (about 20 km s"') with a 
position angle around 80°. The direction and the velocities sug- 
gest that they could be part of the same outflow, generated by the 
same deflection mechanism. 



Groups H, I and J, on the northeast side, are moving NE 
(P.A.~45°) with a Vtan between 10 and 50kms"'. 

Group N, in the central region of the main body, could be part 
of the HH 54 flow, with a vw -20-25 km s"' and a P A.~320°. 

It is difficult to determine the membership of several knots 
(as B or F groups) in the central part of the HH object. The 
bulk of knot B possibly moves slightly towards ESE, while Bl, 
few arcseconds north, seems to proceed in the opposite direction. 
Also F and Fl show a similar behaviour. 

The molecular component of knot B presents distinct inter- 
nal motions. We detect three components, reported in Tab. IB. 41 
(Bl, B2, B3, from north to south). Bi appears as a bright arc 
in the high spatial resolution ISAAC images and it is moving 
towards NNE (~20°). B2 and B3 are proceeding towards ESE 
(130°- 140°) as in the optical. In this case we are possibly ob- 
serving the motion of the two different outflows. 

Appendix C: Flux and velocity variability 

C.1. HH52 

C.1.1. HH52BandC 

Both wings of HH 52 bow-shock show a high degree of variabil- 
ity. As an example of such a variability, we show in Fig. IC.ll 
a close up of HH52 in 1993 and 2006 Ha calibrated images. 
Changes in morphology and flux are easy to recognise, espe- 
cially for the knots of group B. 



A.Caratti o Garatti et al. 



Time-variability in the Clia II outflows, Online Material p 5 



HH 54 




Fig. B.l. Flow charts with error bars of HH 54 (left panel), HH 52 and 53 (right panel) in H2 filter. Proper motions and their error 
bars are indicated by arrows and ellipses, respectively. 




Fig. C.l. Variability in morphology and flux of HH 52. Left panel: close up of HH 52 from 1993 Ha calibrated image (labels refer 
to 2006 image). Right panel: close up of HH52 from 2006 Ha calibrated image. Both images have the same contour levels. 



HH52 Bl and B2 come from the fragmentation of a sin- 
gle structure around 1987. The first knot is moving ESE with a 
P.A.~110°, while the second is moving northward (P.A.~5°). A 
third knot B3, clearly visible from 1989, likely causes the break, 
proceeding from SE to NE and pushing forward the structure. 
The direction of the feature considerably changes with time in 
both optical filters, accordingly, from ~55° in 1989 to 85° in 
1993. Also Bl and B2 show variability (see Fig.|CT] Fig.lH and 
Fig. IC.2b . exhibiting evidence of deceleration and acceleration, 
respectively. 

Knots HH52B4-B9, in the outskirts of the wing, exhibit 
large variability and complicated motions as well. The presence 
of such a high variability and the large gap between the sec- 
ond to last (1993 for [S 11] and 1995 for Ha) and the last epoch 
image (2006) produce uncertainties in the identification of the 
knots and in the P.M. derivation. A large structure, B5, appears 



in 1989 (see Fig.|7] top panel and also Fig. IC.lt , rapidly increas- 
ing its brightness and then almost vanishing in 2006, apparently 
fragmented in several knots (B4, B5, B6, and B8), that move in 
different directions. The flux variation of the object has different 
timing in the three filters, with a delay of some years between 
the flux variations of [S 11] and Ha, and then between Ha and 
H2 (Fig |7] top panel). In [S 11] it slightly varies its brightness, 
reaching a maximum possibly around 1991-1993. In Ha and H2 
the variations are more considerable, especially in the molecu- 
lar component, where the flux increases an order of magnitude 
during six years. 

Knots C in the left wing of HH 52 present a peculiar be- 
haviour as well. CI (see Fig.|2l central panel and Fig. |C.3l l shows 
a large variability similar to B5, with a large brightness increase 
in Ha, due to the presence of a fast shock colliding against 
the target. This is clearly visible in 1993 and 1995 Ha images 
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Fig.C.2. Variability in HH52B2. Bottom panel shows the 
P.M.s as a function of time in [S ii] (circles)and Ha (triangles). 
The dashed line is the best fit to the data points. The slope of 
the fit gives the acceleration of the knot. In the top panel of the 
figure the measured fluxes (uncorrected for the extinction) in the 
two filters are reported. 



Fig.C.4. Variability in HH54Y2. Bottom panel shows the 
P.M.s as a function of time in [S ii] (circles). Ha (triangles), and 
H2 (squares). The dashed line is the best fit to the data points. 
The slope of the fit gives the acceleration of the knot. In the top 
panel of the figure the measured fluxes (uncorrected for the ex- 
tinction) in the three filters are reported. 
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Fig. C.3. Variability in HH52C1. Bottom panel shows the 
P.M.s as a function of time in [S 11] (circles)and Ho' (triangles). 
The dashed line is the best fit to the data points. The slope of 
the fit gives the acceleration of the knot. In the top panel of the 
figure the measured fluxes (uncorrected for the extinction) in the 
two filters are reported. 



(Fig. IC.lt . where a high velocity knot (C6) appears behind CI . In 
the 2006 image the bullet is not visible anymore, but an increase 
in CI brightness is detected, in both atomic emissions. As in B5, 
also the molecular emission increases with time. Here the incre- 
ment is not so pronounced as in the previous knot. HH 52 C2, 
located behind HH52 nucleus, starts forming in 1989 as a [S 11] 
emission (Fig|71 bottom panel) and increases its luminosity until 
2006. On the other hand the Ha emission is visible above a 3cr 
limit only in the 2006 image. 

C.2. HH54 

C.2.1. HH 54 streamer 

Figure IC.4I (bottom panel) displays the peculiar kinematics of 
HH54Y2, located in the middle of the streamer. The proper 
motion measurements reveal an apparent acceleration of the 
knot, detected in all three filters. However, due to the faint 
emission, errors in H2 appear remarkably larger Here, the 
P.M. measurements are referred to epoch 2002, since in H2 
2005 image Y2 is out of the FoV We measure an accelera- 
tion of 0.006+0.0002" yr-2 and 0.01 +0.004" yr ^ in Ha and 
H2, respectively (see Table IC.ll ). The average value obtained 
including both optical (Ha and [Sii]) and molecular P.M.s is 
0.0031 ±0.0008" yr"^. Thus a mechanism must exist there that 
is accelerating the streamer or, conversely, the observed motion 
could be a reminiscence of the original acceleration experienced 
by the flow. The upper panel of Fig. |C. 41 reports the measured 
flux for the three filters. For this object, all the emissions appear 
almost constant with time (inside the error bars). This is indeed 
due to the high degree of ionisation along the HH 54 streamer. 
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Fig. C.5. Measured fluxes (uncorrected for the extinction) in 
knots HH54A1 (bottom panel), and A8 (top panel) in [Sii] 
(circles). Ha (triangles), and H2 (squares) filters. 



Also along the streamer it is possible to observe the forma- 
tion of new condensations. A new Wa knot is forming ~1".6 
ENE of HH54Y3. This is probably an instability in the flow 
created by the material compressed and pushed ahead of Y3. 

C.2.2. HH54A 

Figure |C3] illustrates flux variability in HH 54 A group. Here, at 
variance with other knots, the measured velocities are constant 
(inside the error bars) in all the three filters. The brightest feature 
(knot Al, bottom panel) doubles its He luminosity, after 1990,. 
Variability is, however, not observed in [S n], which appears con- 
stant during the considered period. The molecular emission, that 
does not spatially coincide with the atomic one but it is located 
about 1".6 forward, from 1993 shows an increase with time. 
Knot A8, next to Al along the flow, appears in Hff images in 
1989 (Fig. IC.5l top panel), and increases its brightness in all the 
three filters, and it is extremely luminous in the molecular com- 
ponent. In both [S 11] and H2 it has a diffuse emission, probably 
originating from moving gas cooling behind the Ha emission. 
The feature has a different direction (P.A.~35°) with respect to 
the bulk of the flow (P.A.~80°). This deflection could be gener- 
ated in the interaction with knot A6 (see also Fig. IA.2b . a faint 
emission following the group, that seems to push aside A8, or, 
conversely, it could indicate that the knot is not part of the same 
flow. 

C.2.3. HH54G 

Few arcseconds ahead of knots G and GO (see Section |4!4l also 
knots HH54 Gl and G3 are deflected (P.A.~60-80°) and ac- 
celerated (see Fig. IC.6I IC.7l i. likely by the same mechanism of 
HH54G/G0. 



Fig. C.6. Variability in HH54G1. Bottom panel shows the 
P.M.s as a function of time in [S 11] (circles)and Ha (triangles). 
The dashed line is the best fit to the data points. The slope of 
the fit gives the acceleration of the knot. In the top panel of the 
figure the measured fluxes (uncorrected for the extinction) in the 
two filters are reported. 
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Fig.C.7. Variability in HH54G3. Bottom panel shows the 
P.M.s as a function of time in [S 11] (circles)and Ha (triangles). 
The dashed line is the best fit to the data points. The slope of 
the fit gives the acceleration of the knot. In the top panel of the 
figure the measured fluxes (uncorrected for the extinction) in the 
two filters are reported. 
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deceleration of 0.02+0.01" yr^ in Ha (see Table ICTT i. 
and a value of 0.0025±0.0012"yr^ or 0.0063+0.0005" yr^^ 
without considering 1993 data, combining all the opti- 
cal P.M.s (see upper Fig. IC.8I bottom panel). During 
this period the Ha and [S ii] fluxes gradually increase 
from 4.9x10"'^ ergs"' cm"^ to 7.2xl0"'^ergs"' cm"^, and 
1.8x10"'^ erg S-' cm-2 to 2.8x10"'^ ergs 'cm ^, respectively 
(top panel). 

Knot H2 (see lower Fig. |C.8l bottom panel) apparently drops 
in velocity between 1987 and 1989 Ha, changing the P.M. value 
from -0.08 to 0.04 "yr"' (decelerating of 0.022+0.002" yr^^ 
in Ha, see Table IC.lt . The knot has also a relevant flux 
variabiUty at least in the Ha and H2 emissions, that have 
a constant increase of the flux from 2.7xl0~'^ergs"' cm"^ 

(2006) in Ha, and from 
-2 (2005) 

in H2. On the other hand the measured flux in the three [Sii] 
images is constant. 

Finally H3 emission appears only in 2006 optical images, 
whereas in H2 it appears as a faint featureless emission and no 
significant motion is detectable. Again, the gap between the sec- 
ond to last and the last optical images makes it hard to derive the 
initial position of the condensation. We have tried to correlate it 
with a feature close to H2 in 1993, but the interpretation is not 
certain. 
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Fig. C.8. Variability in HH 54 HI upper Figure and H2 bottom 
Figure The bottom panels show the P.M.s as a function of time 
in both [S 11] (circles) and Ha (triangles). The continuous and 
dashed lines in the top figure are the best fits of the data points 
with and without 1993 data points (see text), respectively. The 
slope of the fits gives the deceleration of the knot. In the top 
panel of the figures the measured fluxes (uncorrected for the 
extinction) in the two filters are reported. 



C.2.4. HH54H 

Group HH 54 H is a very active and fast variable region of the 
flow. 

Between 1987 and 1993 HI has decreased the speed 
in the optical bands from 40 to 20kms-' with a P.M. 



Appendix D: Radial velocity 

D.1. HH 52 flow 

Along the HH 52 flow, knots have radial velocities from about 
-70 to -90km s ', with lower values detected at the end of the 
flow, in HH54, where the velocity decreases moving from G 
to C, indicating that the flow has been decelerated. The lowest 
value is detected in HH54C2 (-18+4kms-') at the very end of 
the flow. C3 is following C2 at a higher velocity (-45+4 km s"') 
(see Tab. ID. 3b . Along the jet, in knots G1-G3, we observe two 
distinct velocity components of -87 and -21 km s ' . In the EMMl 
spectral image the emissions are well separated, delineating two 
distinct lines of similar intensities for each detected species (see 
Tab. ID. 3b . At the base of the jet, where GO and G are located, 
the fast component, extending few arcseconds from the HH52 
streamer to GO, shows an abrupt deceleration (from -78 to - 
18kms ') and then ahead, along the flow, a second slower ve- 
locity component is detected. Also in consideration of the proper 
motion analysis, where a portion of the flow appears to be de- 
flected in this region, the two velocities could be ascribed to 
the two flow components, rather than a bow shock structure. 
Moreover, due to the profile of the observed emission lines, a 
hypothetical bow shock would lie close the plane of the sky (see 
e. g. Hartigan et al. 1 1987b . but, as the inclination angle suggests 
(see also Sect. 14. 6b . this is not the case. 

Two velocities components, not completely separated, are 
also detected in HH52 knots located in front (D4) and 
behind (A3-A4) the nucleus (Al, v™^/=-79±6kms"') (see 
Tab. ID. lb . Velocities in A3-A4 are slightly higher (-98+5 km s"' 
and -35+6 km s"',) with respect to D4 (-85+3 km s"' and - 
30+3 km s"'). Line fluxes of the low velocity component are 
considerably fainter than the other component. The velocity pro- 
files of the Ha line in both structures closely match the theoreti- 
cal profile of a bow shock with an inclination angle of ~60° with 
respect to the plane of the sky (see Hartigan et al. 11987b and 
a shock velocity of 70-80 km s"' (observed from the full width 
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zero intensity - FWZI, measured on the spectrum where the flux 
reaches a 2cr background noise level, see e. g. Davis et al. |2001| l. 

D.2. HH 53 flow 

The highest Vrad values are observed in the three knots of HH 53 
(see Tab.ESl outflow, A (-llO+lOkms-'), B (-113+9 kms^i), 
and C (-96+ 14km s"'), derived averaging the radial velocities 
of the different species from B&C slit 1 and 2. In the EMMI 
spectrum (Tab. ID. 21 ) it is not possible to spatially disentan- 
gle the emission coming from knots C and CI. Here, how- 
ever, we measure two different velocities of -104+7 km s"' and 
-55+18 km s ' . Considering that the two knots are from two sep- 
arate flows, the higher velocity should be ascribed to knot C and 
the lower to C 1 . 

D.3. HH 54 flow 

Along the portion of the HH 54 streamer encompassed by our 
B&C slits, we measure an increase in the velocity moving to- 
wards the main body of HH 54 (see also Tab. |D.3t . The val- 
ues rise from about -50kms"', roughly at the base of the 
streamer (X1A-X3), up to -lOOkms"', close to the middle 
(X4A-X4B). Such a behaviour was also observed by Graham 
& Hartigan ( I1988D . but the value measured in the middle (con- 
sidering also the errors) is lower than ours of ~10-20kms 
Such a difference in their values could be due to a not perfect 
positioning of the slit, due to the sinuous geometry of the HH 54 
streamer. 

Appendix E: Inclination and spatial velocity 

In HH 52 bow shock we obtain an average inclination with re- 
spect to the sky plane of 58°+3°, excluding D2-D3 measure- 
ment, that shows a value of 67°+4°. HH 53 El and F2, as already 
deduced from the P.M. analysis, are part of the same outflow 
with inclinations of 57°+3° and 63°+4°, respectively. The asso- 
ciation of groups HH 54 G and C to the HH 52 outflow is also 
confirmed. HH 54 GO has an inclination of 55°+3°, at the end of 
the streamer Moreover, if we assume that along HH54G1-G3 
the high and low radial velocities are associated with the high 
and low tangential velocities observed along the jet, we obtain 
inclinations ranging between 50° +4° and 56°+5°, respectively. 
Values in HH54 C2 and C3 are 61°+3° and 52°+8°, respec- 
tively. Spatial velocity along the flow are around 100 km s"' (see 
Tab. IE. lb and decreases in group G and C (20-50kms '). 

HH53 A, B, and C have similar inclinations (83°+2°) and 
spatial velocities (around llOkms"'). HH53C1 is superim- 
posed on the second flow and is moving accordingly with the 
HH 52 streamer, with an inclination angle of ~49°+9°, consider- 
ing the low velocity component measured in HH53 C-Cl spec- 
trum (see Tab. ID. 21 ). 

In the HH54 streamer we measure an average inclination 
angle of 67+3°, but the errors of the single data points are too 
large to derive variations of inchnation along the wiggling jet. 
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Table B.l. Proper Motions in HH 52 & HH 53 - Ho- & [S n] - Parameters derived from the linear fit to be continued 
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12 


55 


09.5 


-76 


57 


52 


0.035±0.003 


78±4 


29±3 


Ha 


HH52B1 


12 


55 


05.6 


-76 


57 


54 


0.046±0.002° 


116±3 


39±2° 


[Sii] 




12 


55 


05.2 


-76 


57 


55 


0.031 ±0.004° 


111±6 


26±3° 


Ha 


HH52B2 


12 


55 


05.3 


-76 


57 


53 


0.053±0.002'' 


2±3 


45±2° 


[Sii] 




12 


55 


05.2 


-76 


57 


53 


0.062±0.002'' 


1±3 


53±2° 


Ha 


HH52B3 


12 


55 


04.6 


-76 


57 


54 


0.038±0.002° 


71±3 


32±2° 


[Sii] 




12 


55 


04.6 


-76 


57 


54 


0.044±0.003 


80±3 


37±2 


Ha 


HH52B4 


12 


55 


03.0 


-76 


57 


57 


0.017±0.002'' 


63±7 


14±2° 


[Sii] 




12 


55 


03.0 


-76 


57 


58 


- 


- 




Ha 


HH52B5 


12 


55 


02.9 


-76 


57 


58 


0.011±0.002° 


309±10 


10±2° 


[Sn] 




12 


55 


02.7 


-76 


57 


58 


0.097±0.003 


195±2 


82±2 


Ha 


HH52B6 


12 


55 


02.2 


-76 


57 


55 


0.034±0.002 


257±3 


29±1 


[Sn] 




12 


55 


02.1 


-76 


57 


55 


0.033±0.001 


249±2 


28±1 


Ha 


HH52B7 


12 


55 


03.4 


-76 


58 


00 


0.033±0.003 


157±6 


28±3 


[Sn] 




12 


55 


03.2 


-76 


58 


00 


0.062±0.003 


152±3 


52±2 


Ha 


HH52B8 


12 


55 


02.4 


-76 


57 


57 


0.019±0.005'' 


288±11 


16±4° 


[Sn] 




12 


55 


02.5 


-76 


57 


56 


- 


- 




Ha 


HH52B9 


12 


55 


03.6 


-76 


57 


57 


0.045±0.011 


82±15 


38±9 


Ha 


HH52C1 


12 


55 


05.4 


-76 


58 


01 


0.037±0.002'' 


229±3 


31±1° 


[Sn] 




12 


55 


05.4 


-76 


58 


01 


0.009±0.002 


47±14 


8±2 


Ha 


HH52C2 


12 


55 


05.5 


-76 


57 


58 


0.056+0.002° 


38±2 


48±2° 


[Sn] 




12 


55 


05.4 


-76 


57 


58 


0.046±0.004 


56±5 


39±3 


Ha 


HH52C3 


12 


55 


06.2 


-76 


57 


59 


- 


- 


- 


[Sii],Ha 


HH52C4 


12 


55 


06.7 


-76 


57 


58 


0.056+0.002° 


4±2 


47±2° 


[Sn] 




12 


55 


06.7 


-76 


57 


58 


0.068±0.003° 


4±3 


58±2° 


Ha 


HH52C5 


12 


55 


08.0 


-76 


58 


01 


- 


- 


- 


[Sn] 


HH52C6 


12 


55 


05.3 


-76 


58 


01 


- 


- 


- 


[Sii], Ha 


HH52D1 


12 


55 


13.9 


-76 


57 


47 


- 


- 


- 


[Sn] 




12 


55 


13.9 


-76 


57 


47 


0.060±0.006 


35±5 


50±5 


Ha 


HH52D2 


12 


55 


12.8 


-76 


57 


49 


0.046±0.002 


32±2 


38±2 


[Sn] 




12 


55 


12.7 


-76 


57 


49 


0.049±0.002 


23±3 


41±2 


Ha 


HH52D3 


12 


55 


12.1 


-76 


57 


49 


0.046±0.002 


49±2 


39±2 


[Sn] 




12 


55 


12.0 


-76 


57 


48 


0.047±0.003 


22±3 


40±2 


Ha 


HH 52 D4 


12 


55 


10.0 


-76 


57 


46 


0.044±0.002° 


47±2 


37±1° 


[Sn] 




12 


55 


09.7 


-76 


57 


47 


0.043±0.003 


57±4 


36±3 


Ha 


HH 52 D5 


12 


55 


06.7 


-76 


57 


43 


0.036±0.002° 


53±3 


30±2° 


[Sn] 




12 


55 


06.6 


-76 


57 


44 


0.040±0.003 


24±4 


34±3 


Ha 


HH53 A 


12 


55 


15.4 


-76 


57 


30 


0.019+0.002° 


303±5 


16±1° 


[Sn] 




12 


55 


15.4 


-76 


57 


30 


\J.\jlV+\J.Wi 


Z / J±0 


to±2 


rlQr 


HH53B 


12 


55 


14.5 


-76 


57 


31 


0.007±0.002° 


235±14 


6±2° 


[Sn] 




12 


55 


14.5 


-76 


57 


31 


0.013±0.003 


252±11 


11±2 


Ha 


HH53B1 


12 


55 


13.6 


-76 


57 


31 


0.081 ±0.002 


262±1 


69±2 


Ha 


HH53B2 


12 


55 


14.1 


-76 


57 


31 








[Sii],Ha 


HH53C 


12 


55 


16.5 


-76 


57 


30 


0.013±0.002 


239±9 


11±2 


[Sn] 




12 


55 


16.5 


-76 


57 


30 


0.014±0.003 


232±12 


12±3 


Ha 


HH53C1 


12 


55 


16.1 


-76 


57 


30 


0.060±0.006 


58±6 


50±5 


Ha 



Notes: 

" Mean value, variable in velocity, see text and Table lC.l1 
* Below 3cr. 
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Table B.l. Proper Motions in HH 52 & HH 53 - Ho- & [S n] - Parameters derived from the linear fit continued 



Knot ID 




t!'(2000.0) 


^(2000.0) 


P.M. 


PA. 


^ utn 


Filter 




('• 
\ 


/)( 




(° 




") 




\ } 


(km s"') 




HH53E 


12 


55 


19.4 


-76 


57 


21 


0.016±0.002 


240±6 


14±2 


[Sn] 




12 


55 


19.3 


-76 


57 


21 


- 


~ 


~ 


Ho- 


HH53E1 


12 


55 


19.9 


-76 


57 


24 


0.052±0.004 


63±4 


44±3 


ls n] 




12 


55 


19.8 


-76 


57 


24 


0.049±0.003 


97±3 


41±3 


Ho- 


HH53F 


12 


55 


23.0 


-76 


57 


23 


0.031 ±0.002 


64±3 


26±1 


ls n] 




12 


55 


23.0 


-76 


57 


23 


0.049±0.003 


66±3 


41±3 


Ho 


HH53F1 


12 


55 


24.6 


-76 


57 


20 


0.079+0.002° 


64±2 


67±2'' 


[Sn] 






55 


24.5 


-76 


57 


20 


U.U4y+U.UUo 


oy+i 


4z±2 


rlo 


HH53F2 


12 


55 


22.1 


-76 


57 


28 


0.051 ±0.008 


52±8 


43±6 


Ho 


HH53G 


12 


55 


30.3 


-76 


57 


11 


0.063±0.002'' 


67±2 


53±2" 


[Sn] 




12 


55 


30.3 


-76 


57 


11 


0.037±0.003° 


61±4 


31 ±2° 


Ho 


HH53H 


12 


55 


15.4 


-76 


57 


25 








[Sn] 




12 


55 


15.3 


-76 


57 


25 


0.065±0.004 


52±3 


55±3 


Ho 


HH53I 


12 


55 


13.1 


-76 


57 


34 


0.063±0.002 


51±2 


53±2 


[Sn] 




12 


55 


12.2 


-76 


57 


36 


0.066±0.003 


53±2 


56±2 


Ho 



Notes: 

" Mean value, variable in velocity, see text and Table lC.l1 
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Table B.2. Proper Motions in HH 54 - Ha & [S ii] - Parameters derived from the linear fit to be continued 



Knot ID 




Q'(2000.0) 


(5(2000.0) 


P.M. 


P.A. 


V, 
^ !ati 


Filter 




('' 


m 




(° 




") 






(km s"') 




HH54A1 


12 


55 


49.4 


-76 


56 


30 


0.036+0.004 


70+6 


31+3 


[Sn] 




12 


55 


49.4 


-76 


56 


30 


0.035+0.002 


85+2 


29+2 


Ha 


HH54A2 


12 


55 


50.2 


-76 


56 


30 


0.036+0.004 


81+6 


30+3 


[Sn] 




12 


55 


50.2 


-76 


56 


30 


0.034+0.002 


73+3 


29+2 


Hot 


HH54A3 


12 


55 


48.7 


-76 


56 


33 


0.013+0.004 


166+17 


11+3 


[Sn] 




12 


55 


48.7 


-76 


56 


33 


0.015+0.002 


121+6 


13+3 


Ho- 


HH54A4 


12 


55 


47.7 


-76 


56 


33 


0.030+0.004 


225+7 


25+3 


[Sn] 




12 


55 


47.8 


-76 


56 


34 




- 


- 


Hff 


HH54A5 


12 


55 


47.4 


-76 


56 


30 


0.066+0.006 


125+6 


56+5 


[Sn] 




12 


55 


AHA 


- /o 


56 


31 


- 




- 


Ho- 


HH54A6 


12 


55 


AO ^ 

48.1 


- /o 


56 


1 1 

31 


0.031+0.007 


53+13 


27+6 


[Sn] 




12 


55 


48. i 


-lb 


56 


31 


0.025+0.005 


40+11 


21+4 


Ha 


HH54A7 


12 


55 


49.4 


-76 


56 


32 


0.017+0.005 


88+27 


14+4 


[Sn] 




12 


55 


49.4 


-76 


56 


32 


0.016+0.002 


23+6 


14+1 


Ha 


HH54A8 


12 


55 


48.7 


-76 


56 


29 


0.039+0.009 


28+20 


33+8 


[Sn] 




12 


55 


48.7 


-76 


56 


29 


0.034+0.002 


66+3 


29+2 


Ha 


HH 54 B 


12 


55 


50.9 


-76 


56 


21 


0.008+0.004* 


110+30 


6+3* 


[S n] 




12 


55 


50.9 


-76 


56 


21 


- 


~ 


- 


Ha 


HH54B1 


12 


55 


50.6 


-76 


56 


19 


0.010+0.005* 


282+32 


8+4* 


[Sn] 




12 


55 


50.6 


-76 


56 


19 


- 






Ha 


HH54B2 


12 


55 


49.8 


-76 


56 


21 




- 




[Sn] 




12 


55 


49.8 


-76 


56 


22 


0.017+0.002 


313+7 


15+2 


Ha 


HH54B3 


12 


55 


50.5 


-76 


56 


24 




- 


- 


Ha 


HH54C1 


12 


55 


51.3 


-76 


56 


05 


0.032+0.004 


73+8 


27+3 


[S n] 




12 


55 


51.3 


-76 


56 


05 


0.017+0.002 


67+5 


14+1 


Ha 


HH54 C2 


12 


55 


50.4 


-76 


56 


03 


0.041+0.004 


42+5 


34+3 


[Sn] 




12 


55 


50.4 


-76 


56 


03 


0.017+0.002 


70+5 


14+1 


Ha 


HH54C3 


12 


55 


51.0 


-76 


56 


02 


0.017+0.006* 


239+19 


14+5* 


[Sn] 




12 


55 


51.0 


-76 


56 


02 


0.016+0.002" 


267+6 


13+2° 


Ha 


HH54E 


12 


55 


53.7 


-76 


56 


22 


0.017+0.005 


35+15 


14+4 


[Sn] 




12 


55 


53.7 


-76 


56 


22 


0.018+0.002 


16+6 


15+2 


Ha 


HH54F 


12 


55 


49.0 


-76 


56 


19 


0.018+0.003 


125+11 


15+3 


[Sn] 




12 


55 


49.0 


-76 


56 


19 


0.013+0.002 


129+7 


11+2 


Ha 


HH 54F1 


12 


55 


48.4 


-76 


56 


17 


0.015+0.002" 


251+9 


13+2" 


Ha 


HH 54 G 


12 


55 


45.4 


-76 


56 


13 


0.048+0.004" 


78+4 


41+3" 


[S n] 




12 


55 


45.4 


-76 


56 


13 


0.035+0.001" 


71+3 


29+1" 


Ha 


HH54 GO 


12 


55 


44.8 


-76 


56 


13 


0.036+0.005" 


36+7 


31+4" 


[Sn] 




12 


55 


44.8 


-76 


56 


13 


0.055+0.003" 


52+3 


46+3" 


Ha 


HH54G1 


12 


55 


47.3 


-76 


56 


09 


0.096+0.006" 


81+3 


81+5" 


[Sn] 




12 


55 


47.3 


-76 


56 


09 


0.073+0.006" 


72+6 


62+5" 


Ha 


HH54 G2 


12 


55 


47.6 


-76 


56 


05 


0.017+0.005* 


58+17 


14+5* 


[Sn] 




12 


55 


47.8 


-76 


56 


05 


0.079+0.003" 


68+2 


67+2" 


Ha 


HH54 G3 


12 


55 


49.1 


-76 


56 


06 


0.079+0.006" 


63+4 


67+5" 


[Sn] 




12 


55 


49.3 


-76 


56 


06 


0.080+0.002" 


66+2 


68+2" 


Ha 


HH54H1 


12 


55 


52.8 


-76 


56 


05 


0.026+0.004 


36+8 


22+3 


[S n] 




12 


55 


52.8 


-76 


56 


05 


0.033+0.002" 


53+4 


28+2" 


Ha 


HH54 H2 


12 


55 


53.1 


-76 


56 


07 


0.035+0.004 


23+6 


29+3 


[Sn] 




1 1 


55 


53.1 


-76 


56 


07 


U.Uji+U.UUZ 




4o±Z 


rlQr 


HH54 H3 


12 


55 


52.4 


-76 


56 


10 


0.025+0.003 


187+9 


21+3 


[S n] 




12 


55 


52.6 


-76 


56 


11 


0.043+0.003" 


127+4 


37+3" 


Ha 


HH54I 


12 


55 


54.4 


-76 


56 


08 


0.031+0.006 


48+11 


26+5 


[S n] 




12 


55 


54.5 


-76 


56 


08 


0.031+0.002 


41+4 


26+2 


Ha 


HH54J 


12 


55 


52.3 


-76 


56 


17 


0.012+0.005* 


45+25 


10+4* 


[S n] 




12 


55 


52.3 


-76 


56 


18 








[S n],Ha 


HH 54 Jl 


12 


55 


51.4 


-76 


56 


16 








[S n],Ha 



Notes: 

" Mean value, variable in velocity, see text and Table lC.ll 
* Below 3cr. 
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Table B.2. Proper Motions in HH 54 - Ha & [S ii] - Parameters derived from the linear fit continued 



Js^not iU 


Q'(2000.0) 


5(2000.0) 


F.M. 


r.A. 




Filter 




(" 


m 




(° 




") 


("yr 1) 


(°) 


(kms ) 




HH54K 


12 


55 


53.8 


-76 


56 


25 


0.017±0.003 


2±14 


14±3 


[Sii] 




12 


55 


53.8 


-76 


56 


25 


0.021 ±0.003 


329±11 


18±3 


Ha 


HH54K1 


12 


55 


54.1 


-76 


56 


27 


0.042±0.004 


207±5 


36±3 


[Sn] 




12 


55 


54.2 


-76 


56 


27 


0.048±0.005 


233±7 


41±4 


Ho- 


HH54K2 


12 


55 


53.3 


-76 


56 


24 


0.022±0.005 


45±14 


18±4 


ls II] 




12 


55 


53.4 


-76 


56 


24 


0.016±0.002 


345±10 


13±1 


Ho- 


HH 54 M 


12 


55 


51.2 


-76 


56 


08 


0.026±0.004 


76±9 


22±4 


ls II] 




12 


55 


51.0 


-76 


56 


09 


0.019±0.002'' 


81±7 


16±2" 


Ho 


HH54M1 


12 


55 


50.2 


-76 


56 


11 


„ 


_ 


_ 


[Sii], Hot 


HH54N 


12 


55 


48.8 


-76 


56 


23 


0.029+0.005 


312±9 


25 ±4 


[Sii] 




12 


55 


49.0 


-76 


56 


23 


0.021±0.002 


326±6 


22±2 


Ho 


HH540 


12 


55 


47.4 


-76 


56 


22 


0.013±0.004 


328±16 


11±3 


[Sii] 




12 


55 


47.3 


-76 


56 


22 


0.051 ±0.005 


58±6 


43±4 


Ho 


HH54 01 


12 


55 


46.7 


-76 


56 


24 


0.046±0.004 


15±7 


39±3 


[Sii] 




12 


55 


46.6 


-76 


56 


24 


0.037±0.003'' 


47±5 


31±3'' 


Ho 


HH54P 


12 


55 


45.2 


-76 


56 


34 


_ 


_ 


_ 


[Sii], Ha 


HH54Q1 


12 


55 


47.2 


-76 


56 


43 


0.044±0.006 


23±8 


37±5 


[Sii] 




12 


55 


52.5 


-76 


56 


43 


0.041±0.002'' 


37±3 


34±2'' 


Ho 


HH54 Q2 


12 


55 


51.8 


-76 


56 


40 


0.046±0.013 


34±17 


39±11 


[Sn] 




12 


55 


51.8 


-76 


56 


40 


0.043±0.003 


79±4 


37±3 


Ho- 


HH54 Q3 


12 


55 


51.3 


-76 


56 


36 


_ 


_ 


_ 


[Sn] 




12 


55 


51.3 


-76 


56 


36 


0.043±0.003'' 


80±4 


37±2'' 


Ho 


HH54R 


12 


55 


56.3 


-76 


57 


20 


_ 


_ 


_ 


[Sii],Ha 


HH54X0 


12 


55 


40.3 


-76 


57 


50 


0.022±0.003 


325±6 


19±2 


Ho 


HH54X1 A 


12 


55 


44.6 


-76 


57 


35 


0.035±0.003" 


101±4 


29±2" 


Ho 


HH54X1B 


12 


55 


42.8 


-76 


57 


36 


0.013±0.002 


313±11 


11±2 


[Sn] 




12 


55 


42.8 


-76 


57 


37 


0.012±0.003 


275±8 


10±3 


Ho 


HH54X1C 


12 


55 


41.6 


-76 


57 


35 


0.036±0.002'' 


340±4 


30±2° 


Hot 


HH54X1D 


12 


55 


43.8 


-76 


57 


27 


0.037±0.004 


33±7 


31±4 


Ha 


HH54X2A 


12 


55 


40.7 


-76 


57 


42 


0.018±0.006 


204±23 


15±5 


[Sn] 




12 


55 


40.9 


-76 


57 


39 


0.021±0.002 


55±6 


18±2 


Ha 


HH54X2B 


12 


55 


41.3 


-76 


57 


39 


0.013±0.002 


314±11 


11±2 


[Sn] 




12 


55 


41.0 


-76 


57 


41 


0.036±0.002 


253±5 


30±2 


Ho 


HH54X3 


12 


55 


39.8 


-76 


57 


36 


_ 


_ 


_ 


[Sn] 




12 


55 


39.8 


-76 


57 


36 


0.021 ±0.003 


256±13 


18±4 


Ho 


HH54X4A 


12 


55 


40.5 


-76 


57 


28 


0.068±0.007 


45 ±6 


57±7 


Ho 


HH54X4B 


12 


55 


42.0 


-76 


57 


22 


0.043±0.003 


22±5 


37±3 


Ha 


HH54X4C 


12 


55 


41.8 


-76 


57 


25 


0.067±0.006 


49±5 


57±5 


[Sn] 




12 


55 


41.8 


-76 


57 


25 


0.047±0.003 


44±3 


39±2 


Ha 


HH54Y1 




j5 


44. y 


- lb 


J / 


i 1 


0.033±0.005 


33±8 


28±4 


[Sn] 




12 


55 


44.9 


-76 


57 


1 1 


0.018±0.002 


356±5 


16±2 


Ho 


HH54 Y2 


12 


55 


44.7 


-76 


57 


06 


0.038±0.005 


18±9 


32±4 


[Sn] 




12 


55 


44.6 


-76 


57 


07 


0.037±0.002" 


26±3 


31 ±2° 


Ho 


HH54 Y3 


1 1 


JJ 


A1 ^ 

4i.Z 


-lb 


J / 


1 1 


0.042±0.004 


63±5 


36±4 


[Sn] 




1 1 


cc 

JJ 


43. Z 


-lb 


J 1 


1 T 


0.068±0.003" 


43±3 


58±3 


Ho 


HH54Y4 


1 o 

tz 


JJ 


4j.U 


- lb 


J / 


UU 


0.043±0.008 


62±17 


36±7 


[Sn] 




12 


55 


45.0 


-lb 


57 


00 


0.067±0.008 


41±7 


57±7 


Ha 


HH54 Y5 


12 


55 


49.1 


-76 


57 


04 


0.013±0.004'' 


14±19 


11±3'' 


[Sn] 




12 


55 


49.0 


-76 


57 


05 


0.022±0.008'' 


90±20 


19±7'' 


Ha 


HH54 Y6 


12 


55 


45.5 


-76 


56 


53 


0.059±0.006 


69±6 


50±5 


[Sn] 




12 


55 


45.5 


-76 


56 


53 


0.059±0.010 


16±18 


50±9 


Ha 


HH54 Y7 


12 


55 


46.1 


-76 


57 


00 


0.072±0.017 


322±15 


61±14 


Ha 


HH54 Y8 


12 


55 


45.3 


-76 


57 


09 


0.026±0.010'' 


22±11 


22±8* 


Ha 


HH 54 Z 


12 


55 


47.2 


-76 


56 


43 


0.037±0.004 


344±7 


31±3 


[Sn] 




12 


55 


47.2 


-76 


56 


45 


0.010±0.002 


325±13 


8±2 


Ha 


HH 54 Zl 


12 


55 


46.1 


-76 


56 


45 


0.013±0.002 


292±9 


11±2 


Ha 



Notes: 

° Mean value, variable in velocity, see text and Table lC.l1 
* Below 3cr. 
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Table B.3. Proper Motions in HH 52 & HH 53 - 



Knot ID 




a(2000.0) 


5(2000.0) 


P.M. 


PA. 




Filter 




('■ 


m 




(° 




") 


("yr->) 


o 


(kms-') 




HH52B2/B3 


12 


55 


05.0 


-76 


57 


54 


0.026+0.004" 


40+9 


22+4" 


Hz 


HH52B5/B6 


12 


55 


02.7 


-76 


57 


57 


0.041+0.014* 


281+12 


35+12* 


Ho 


HH52C1 


12 


55 


05.4 


-76 


58 


01 


0.008+0.004* 


51+26 


7+3* 


Hz 


HH53 A 


12 


55 


15.4 


-76 


57 


30 


0.016+0.003 


298+15 


14+3 


H2 


HH53B 


12 


55 


14.5 


-76 


57 


31 


0.007+0.003* 


231+34 


6+3* 


H2 


HH53C 


12 


55 


16.5 


-76 


57 


30 


0.010+0.003 


270+26 


9+2 


H2 


HH53H 


12 


55 


15.2 


-76 


57 


26 


0.024+0.004 


27+8 


21+4 


H2 



Notes: 

" Mean value, variable in velocity, see text and Table lC.l1 
* Below 3o-. 



Table B.4. Proper Motions in HH 54 - H2 



Knot ID 




ci'(2000.0) 


5(2000.0) 


P.M. 


PA. 




Filter 




(" 


m 




(° 




") 


("yr-') 


n 


(kms-') 




HH54Ai 


12 


55 


47.3 


-76 


56 


29 


0.018+0.005 


343+16 


15+4 


Ht 


HH54A2 


12 


55 


50.0 


-76 


56 


29 


0.011+0.003 


51+17 


9+3 


H2 


HH54A3 


12 


55 


49.2 


-76 


56 


29 


0.018+0.004 


300+10 


16+3 


H2 


HH54A4 


12 


55 


50.9 


-76 


56 


30 


0.014+0.006* 


62+16 


12+5* 


Ht 


HH54B, 


12 


55 


51.0 


-76 


56 


20 


0.016+0.003" 


15+13 


14+2" 


H2 


HH54B2 


12 


55 


51.1 


-76 


56 


21 


0.008+0.003* 


131+23 


6+3* 


H2 


HH54B3 


12 


55 


51.2 


-76 


56 


22 


0.022+0.003" 


140+9 


18+3" 


H2 


HH54C1 


12 


55 


51.9 


-76 


56 


05 


0.015+0.003 


77+9 


13+3 


H2 


HH54C3 


12 


55 


51.3 


-76 


56 


01 


0.009+0.004* 


239+19 


8+3* 


H2 


HH 54 E 


12 


55 


53.9 


-76 


56 


22 


0.012+0.003" 


48+14 


11+3 


H2 


HH54H2 


12 


55 


53.2 


-76 


56 


07 


0.020+0.003 


39+9 


17+3 


H2 


HH 541 


12 


55 


55.2 


-76 


56 


06 


0.016+0.004" 


3+14 


14+3" 


H2 


HH 541 


12 


55 


53.0 


-76 


56 


17 


0.019+0.003 


45+10 


16+3 


H2 


HH54I1 


12 


55 


52.6 


-76 


56 


16 


0.008+0.004* 


32+36 


8+4* 


H, 


HH 54 K 


12 


55 


54.1 


-76 


56 


25 


0.012+0.003 


19+17 


10+2 


H2 


HH54K1 


12 


55 


54.6 


-76 


56 


27 


0.019+0.005 


243+16 


16+4 


H2 


HH54M 


12 


55 


51.5 


-76 


56 


12 


0.020+0.007* 


110+13 


17+6* 


H2 


HH 54 Ml 


12 


55 


50.3 


-76 


56 


13 


0.021+0.003 


211+8 


18+3 


H. 


HH 54 T 


12 


55 


41.6 


-76 


57 


34 








H2 


HH54Y1 


12 


55 


45.0 


-76 


57 


10 


0.016+0.008* 


67+25 


14+7* 


H2 


HH54Y2 


12 


55 


44.4 


-76 


57 


06 


0.020+0.005" 


21+7 


17+4" 


H2 


HH 54 Z 


12 


55 


47.5 


-76 


56 


45 


0.018+0.004 


39+12 


15+3 


H, 


HH54Z1 


12 


55 


46.3 


-76 


56 


45 


0.015+0.004 


317+13 


13+3 


H7 
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Table C.l. Accelerated Motions in HH 52, HH 53, and HH 54 - Hor, [S n] and H2 - Parameters derived from the quadratic fit 



Knot ID 




P.M. 


Accelerated motion 


P.A.(v) 


P.A.(a) 


Filter 


Note 




xlOO (%) 


("yr-i) 


(" vr"^) 


(°) 


D 




Reliable 


HH52A5 


0.10 


0.071 ±0.0 12 


0.0067±0.0013 


110±8 


133±11 


[Sii] 


no 




1.00° 


0.137±0.086* 


0.018±0.008* 


5±33 


340±26 


Ho- 




HH52B1 


0.05 


0.132+0.017 


0.0226±0.0018 


116±8 


299±5 


ls 11] 


yes 




1.00° 


0.093+0.068*" 


0.02±0.01* 


111±20 


285±16 


Ha 




HH52B2 


0.99 


0.099±0.021 


0.0060±0.0021* 


7±11 


12±23 


[Sn] 


yes 




1.00 


0.131±0.023 


0.0109±0.0027 


9±16 


16±23 


Ho- 




HH52B3 


0.34 


0.122±0.017 


0.0121±0.0019 


112±7 


128±9 


es 11] 


no 


HH 52B2/B3 


1.00° 


0.194±0.077* 


0.061±0.015 


357±22 


348±23 


H, 


no 


HH52B4 


0.22 


0.124±0.016 


0.0134±0.0019 


90±6 


94±7 


[Sii] 


no 


HH52B5 


0.05 


0.106±0.041* 


0.0117±0.0030 


243±13 


238±13 


[Sii] 


no 




0.10 


0.146+0.045 


0.0076±0.0052* 


202±12 


215±36 


Ha 




HH52B8 


1.00 


0.164±0.067* 


0.020±0.007* 


290±16 


291±19 


[Sn] 


yes 


HH52C1 


0.05 


0.052±0.012 


0.0105±0.0014 


69±13 


229±10 


[Sii] 


yes 


HH52C4 


0.28 


0.144±0.025 


0.0105±0.0018 


347±7 


338±11 


[Sii] 


yes 




1.00° 


0.166+0.047 


0.0163±0.0050* 


348±14 


338±21 


Ho 




HH 52 D4 


0.18 


0.073+0.014 


0.0038±0.0013* 


33±9 


16±21 


[Sii] 


no 




0.53 


0.093±0.050'' 


0.011±0.0050* 


12±22 


347±28 


Ha 




HH53A 


0.93 


0.034+0.012* 


0.0040±0.0012 


230±18 


197±20 


[Sn] 


no 


HH53B 


0.96 


0.038±0.022* 


0.0047±0.0012 


332±18 


343±18 


[Sn] 


no 


HH53G 


0.07 


0.039±0.013 


0.0049±0.0014 


104±16 


212±18 


[Sn] 


no 




1.00" 


0.212±0.069* 


0.0276±0.0089 


105±15 


113±18 


Ha 




HH 54Bi 


0.34 


0.04+0.01 


0.0088±0.0024 


40±12 


61±15 


H, 


no 


HH 54 B3 


0.58 


0.067±0.013 


0.0145±0.0023 


173±8 


185±10 


H2 


no 


HH 54 C3 


0.85 


0.046±0.020* 


0.054±0.014 


184±21 


165±20 


Ho 


no 


HH54E 


0.86 


0.027±0.008 


0.0039±0.0015* 


28±19 


16±41 


H2 


yes 


HH 54 Fl 


0.09 


0.072±0.024 


0.084±0.019 


312±12 


324±13 


Ho 


no 


HH54G 


1.00° 


0.134+0.057* 


0.013±0.006* 


67±21 


61±33 


[Sn] 


yes 




0.08 


0.155±0.015 


0.0146±0.0017 


85±7 


89±9 


Ha 




HH54G0 


1.00° 


0.2±0.1* 


0.020±0.007* 


28±16 


220±20 


[Sn] 


yes 




0.12 


0.05±0.02«' 


0.012±0.002 


40±20 


230±13 


Ha 




HH54G1 


1.00" 


0.239+0.067 


0.02±0.009* 


79±15 


76±27 


[Sii] 


yes 




0.06 


0.214+0.031 


0.0177±0.0039 


85±12 


85±19 


Ho 




HH54G2 


1.00° 


0.12+0.11* 


0.017±0.0077* 


70±35 


316±26 


[Sii] 


no 




0.37 


0.125+0.019 


0.0097±0.0026 


70±13 


143±19 


Ha 




HH54G3 


1.00" 


0.291 ±0.059 


0.0299±0.0077 


76±11 


80±15 


[Sn] 


yes 




0.10 


0.134±0.025 


0.0068±0.0025'' 


66±16 


67±38 


Ha 




HH54H1 


0.49 


0.033±0.013* 


0.0063±0.0017 


55±30 


199±17 


Ha 


yes 


HH54H2 


0.10 


0.158±0.017 


0.0228±0.0019 


41±5 


256±4 


Ho 


yes 


HH54H3 


1.00° 


0.062±0.023'' 


0.012±0.004 


127±4 


273±16 


Ho 


no 




0.50 


0.079±0.017 


0.007±0.003* 


66±12 


140±35 


H2 




HH54M 


1.00 


0.071±0.019 


0.069±0.021 


70±14 


67±19 


Ho 


yes 


HH54Q1 


1.00° 


0.165±0.103'' 


0.0171±0.0078* 


12±20 


8±28 


[Sii] 


yes 




0.97 


0.093±0.023 


0.0084±0.0023 


7±22 


349±30 


Ha 




HH54Q3 


0.99 


0.097±0.018 


0.0069±0.023 


79±12 


79±22 


Ha 


yes 


HH54X1 A 


0.57 


0.079±0.015 


0.055±0.018 


105±9 


107±16 


Ha 


yes 


HH54X1C 


0.38 


0.025±0.012* 


0.0065±0.0015 


225±30 


186±14 


Ha 


no 


HH54Y2 


1.00 


0.082±0.014 


0.0057±0.0013 


25±7 


23±12 


Ho 


yes 




0.59 


0.071 ±0.026* 


0.0098±0.0037* 


19±16 


40±20 


H2 




HH54Y3 


0.06 


0.145±0.026 


0.0097±0.0026 


37±9 


32±18 


Ha 


yes 



Notes: 

° Only 3 points in the fit. 
* Below 3a-. 
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Table D.l. Radial velocities and line intensities for different optical lines of individual knots in HH52 from EMMI and B&C 
spectra. Radial velocities are corrected for the cloud speed with respect to the LSR (v/,ss=2kms Knee |1992] l. When detected, 
two velocity components are reported. 



HH52 



Knot 


Line 


A 


Vrad 


F± AF 




Id. 


(A) 


(km s-') 


(lO^'^ergs-'cm-^) 


EMMI 


Al 


[Oi] 


6300.3 


-78±4 


8.1±0.1 




[Oi] 


6363.8 


-77±4 


2.9±0.1 




[Nil] 


6548.1 


-78±4 


2.2±0.1 




Ha 


6562.8 


-70±4 


16.7+0.2 




[Nil] 


6583.4 


-79±4 


5.8+0.1 




[Sii] 


6716.4 


-87±4 


12.0+0.1 




[Sii] 


6730.8 


-86±4 


12.4+0.1 


A3-A4 


[Oi] 


6300.3 


-89±4 ; -43+4 


1.5+0.1 ; 1.4+0.1 




[Oi] 


6363.8 




1.0+0.1 




[Nil] 


6548.1 


-96+4 ; -27±4 


1.5+0.1 ; 0.5+0.1 




Hot 


6562.8 


-101±4 ; -40±4 


2.3+0.2 ; 1.9+0.2 




[Nil] 


6583.4 




0.3+0.1 




[Sii] 


6716.4 


-100±4 ; -34±4 


2.1+0.1 ; 1.3+0.1 




[Sii] 


6730.8 


-102±4 ; -31±4 


1.8+0.2 ; 1.0+0.1 


D4 


[Oi] 


6300.3 


-84±4 ; -30±4 


0.7+0.1 ; 0.3+0.1 




[Nil] 


6548.1 


-84+4 ; -33±4 


0.8+0.1 ; 0.4+0.1 




Hot 


6562.8 


-73±4'' 


2.4+0.2 




[Nil] 


6583.4 




0.5+0.2* 




[Sii] 


6716.4 


-87±4 ; -27±4 


1.4+0.1 ; 0.5+0.1 




[Sii] 


6730.8 


-87+4 ; -29±4 


1.1+0.1 ; 0.5+0.1 


B&C 


A--' 


[Oi] 


6300.3 


-102±20 


40.0+0.2 




[Oi] 


6363.8 


-103±20 


14.3+0.2 




[Nil] 


6548.1 


-94±20 


10.3+0.2 




Hot 


6562.8 


-101±20 


95.6+0.3 




[Nil] 


6583.4 


-89±20 


31.9+0.2 




[Sii] 


6716.4 


-81 ±20 


62.8+0.3 




[Sn] 


6730.8 


-78±20 


58.7+0.3 




[Oi] 


6300.3 


-77±20 


5.3+0.4 




[Oi] 


6363.8 


-83±20 


2.2+0.3 




[Nil] 


6548.1 


-85+20 


2.6+0.2 




Ha 


6562.8 


-85±20 


17.2+0.2 




[Nil] 


6583.4 


-85±20 


8.1+0.2 




[Sii] 


6716.4 


-63±20 


6.5+0.2 




[Sii] 


6730.8 


-67±20 


6.8+0.2 


D2-D3 


[Oi] 


6300.3 


-112±20 


7.7+0.2 




[Oi] 


6363.8 


-132±20 


3.5+0.2 




[Nil] 


6548.1 


-112±20 


1.4+0.2 




Hot 


6562.8 


-94±20 


25.3+0.3 




[Nil] 


6583.4 


-103±20 


5.3+0.3 




[Sii] 


6716.4 


-103±20 


13.1+0.3 




[Sii] 


6730.8 


-99±20 


10.6+0.2 



Notes: 

" Not possible to deblend. 
''S/N between 2 and 3. 
' Blend of knots A1,A3,A4,A6,A7. 
Blend of knots B5,B8,B9. 
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Table D.2. Radial velocities and line intensities for different optical lines of individual knots in HH53 from EMMI and B&C 
spectra. Radial velocities are corrected for the cloud speed with respect to the LSR (v/,ss=2kms Knee |1992] l. When detected, 
two velocity components are reported. 



HH53 



Knot 


Line 


A 


Vrarf 


F + AF 




Id. 


(A) 


(kms-') 


(lO^'^ergs"' cm-^) 


EMMI 


C-Cl 


[Oi] 


6300.3 


-101±4 ; -32+4 


1.9+0.1 ; 0.9+0.2 




[Oi] 


6363.8 




1.0+0.4 




Ho- 


6562.8 


-97±4 ; -47±4 


9.8+0.1 ; 2.4+0.2 




[Nil] 


6583.4 


-100±4 ; -51+4 


1.5+0.1 ; 0.4+0.2 




[Sn] 


6716.4 


-111±4 ; -72±4 


2.2+0.1 ; 0.7+0.1 




[Sn] 


6730.8 


-112±4 ; -74±4 


2.0+0.1 ; 0.7+0.1 


El 


Ho- 


6562.8 


-65±4 


0.9+0.1 




[Sn] 


6716.4 


-74+4 


0.7+0.1 




[Sn] 


6730.8 


-76±4 


0.7+0.1 


B&C slit 1 


A 


[Oi] 


6300.3 


-140+30 


4.2+0.8 




Ha 


6562.8 


-114+30 


25.0+0.6 




[Nn] 


6583.4 


-104+30 


5.3+0.8 




[Sn] 


6716.4 


-117+30 


6.6+0.4 




[Sn] 


6730.8 


-128+30 


6.4+0.5 


B 


[Oi] 


6300.3 


-130+30 


6.3+0.6 




[Nn] 


6548.1 


-141+30 


2.5+0.6 




Ha 


6562.8 


-124+30 


34.3+0.6 




[Nn] 


6583.4 


-103+30 


3.6+0.5 




[Sn] 


6716.4 


-139+30 


10.5+0.5 




[Sn] 


6730.8 


-127+30 


9.0+0.6 


C 


Ha 


6562.8 


-102+30 


8.8+0.6 




[Nn] 


6583.4 


-103+30 


2.7+0.5 




[Sn] 


6716.4 


-101+30 


2.1+0.6 


B&C slit 2 


A 


[Oi] 


6300.3 


-96+20 


6.7+0.1 




[Oi] 


6363.8 


-112+20 


2.2+0.1 




[Nn] 


6548.1 


-76+20 


1.6+0.1 




Ha 


6562.8 


-106+20 


26.3+0.1 




[Nn] 


6583.4 


-82+20 


5.1+0.1 




[Sn] 


6716.4 


-115+20 


7.2+0.1 




[Sn] 


6730.8 


-111+20 


7.8+0.1 


B 


[Oi] 


6300.3 


-97+20 


9.7+0.1 




[Oi] 


6363.8 


-96+20 


2.9+0.1 




[Nn] 


6548.1 


-94+20 


2.1+0.1 




Ha 


6562.8 


-98+20 


33.3+0.1 




[Nn] 


6583.4 


-94+20 


6.9+0.1 




[Sn] 


6716.4 


-106+20 


11.5+0.1 




[Sn] 


6730.8 


-112+20 


11.3+0.1 


C 


[Oi] 


6300.3 


-109+20 


3.1+0.1 




[Oi] 


6363.8 


-125+20 


1.1+0.1 




[Nn] 


6548.1 


-120+20 


0.6+0.1 




Ha 


6562.8 


-122+20 


14.2+0.1 




[Nn] 


6583.4 


-120+20 


2.0+0.1 




[Sn] 


6716.4 


-116+20 


4.5+0.1 




[Sn] 


6730.8 


-111+20 


3.9+0.1 


F2 


[Oi] 


6300.3 


-87+20 


1.5+0.1 




[Oi] 


6363.8 




0.8+0.4 




[Nn] 


6548.1 




1.0+0.2 




Ha 


6562.8 


-97+20 


7.3+0.2 




[Nn] 


6583.4 


-85+20 


1.7+0.2 




[Sn] 


6716.4 


-81+20 


2.1+0.2 




[Sn] 


6730.8 


-95+20 


1.7+0.2 



Notes: 

" Not possible to deblend. 
''S/N between 2 and 3. 
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Table D.3. Radial velocities and line intensities for different optical lines of individual knots in HH 54 from EMMI and B&C 
spectra. Radial velocities are corrected for the cloud speed with respect to the LSR (v/,ss=2kms Knee |1992] l. When detected, 
two velocity components are reported. 

HH54 



Knot 


Line 


A 




F±AF 




Id. 


(A) 


(km s-') 


(I0"'^ergs-'cm-2) 


EMMI 



C2 


[Oi] 


6300.3 


-40±4 


6.4+0.3 




[Oi] 


6363.8 


-43±4 


2.8+0.3 




[Nil] 


6548.1 


-47±4 


2.2+0.2 




Ha 


6562.8 


-52±4 


13.9+0.3 




[Nil] 


6583.4 


-44±4 


6.5+0.2 




[Sii] 


6716.4 


-46±4 


8.9+0.2 




[Sii] 


6730.8 


-46±4 


8.9+0.2 


C3 


[Oi] 


6300.3 


-12±4 


1.5+0.1 




[Oi] 


6363.8 


-17±4 


0.7+0.2 




[Nil] 


6548.1 


-16+4 


0.4+0.1 




Ha 


6562.8 


-19±4 


2.8+0.1 




[Nil] 


6583.4 


-20±4 


1.6+0.1 




[Sii] 


6716.4 


-17±4 


1.7+0.1 




[Sii] 


6730.8 


-22±4 


1.8+0.1 


G-GO 


[Oi] 


6300.3 


-71+4 ; -24+4 


0.4+0.1 ; 0.6+0.1 




[Nil] 


6548.1 




0.5+0.2 ; - 




Ha 


6562.8 


-80±4 ; -15±4 


3.0+0.1 ; 0.5+0.1 




[Nil] 


6583.4 


-73±4 ;-15±4 


0.8+0.1 ; 0.6+0.1 




[Sii] 


67164 


-78±4 ; -17±4 


1.0+0.2 ; 0.6+0.2 




[Sii] 


6730.8 


-79±4 ;-21±4 


1.0+0.2 ; 0.5+0.2 


G1-G3" 


[Oi] 


6300.3 


-86±4 ;-18±4 


1.9+0.1 ; 2.5+0.1 




[Oi] 


6363.8 




0.6+0.2 ; 0.8+0.2 




[Nil] 


6548.1 




0.6+0.2 ; 1.0+0.3 




Ha 


6562.8 


-82±4 ; -24±4 


5.6+0.3 ; 5.0+0.3 




[Nil] 


6583.4 


-86±4 ; -20±4 


1.2+0.2 ; 2.1+0.2 




[Sii] 


67164 


-89±4 ; -20±4 


2.6+0.2 ; 2.3+0.2 




[Sii] 


6730.8 


-91±4 ;-22±4 


1.6+0.2 ; 2.0+0.1 


B&C slit 3 


XID-X4C 


[Oi] 


6300.3 


-34±30 


6.5+1.3 




Ha 


6562.8 


-70±30 


40.5+2.8 




[Nil] 


6583.4 


-58±30 


5.4+2.0 




[Sii] 


6716.4 


-61±30 


6.6+0.6 




[Sii] 


6730.8 


-67±30 


6.3+2.5 


XIA 


Ha 


6562.8 


-70±20 


4.5+0.1 


X3 


[Oi] 


6300.3 


-40±20 


1.6+0.1 




Ha 


6562.8 


-42±20 


9.9+0.1 




[Nil] 


6583.4 


-65±20 


1.1+0.1 




[Sii] 


6716.4 


-53±20 


2.0+0.2 




[Sii] 


6730.8 


-34±20 


1.3+0.2 


X4A-X4B 


Ha 


6562.8 


-115+20 


11.8+0.2 




[Nil] 


6583.4 


-110+20 


0.8+0.2 




[Sii] 


6716.4 


-92+20 


1.8+0.2 




[Sii] 


6730.8 


-82+20 


1.9+0.2 



Notes: 

" Blend of knots G1,G2 and G3. 
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Table D.4. [Fe n] radial velocities of the individual knots of HH 54 obtained from ISAAC high resolution spectroscopy, corrected 
for the cloud speed with respect to the LSR (v/,5R=2kins Knee |1992l l. When detected, two velocities components are reported. 



HH54 


- ISAAC slit 1 


knot 


v„rf (kms"') 


Al 


-52±3 ; -105±3 


A3 


-34±3 


B 


-14±3 ; -109±3 


B3 


-49±3 ;-105±3 


J-Jl 


-9±3 ; -91±3 


H2 


-36±3 ; -90±3 


H3 


-20±3 ; -45±3 


Z 


-100+3 


HH54 


- ISAAC slit 2 


knot 


v,ad (kms"') 


CI 


-5±3 ; -90±3 


C3 


-9±3 ; -85±3 


E 


-52+3 


J-Jl 


-18±3 ;-62±3 


H3 


-40±3 


K 


-62±3 


Kl 


-45±3 



Table E.l. Inclination and spatial velocities of individual knots of HH 52, 53, and 54 inferred from the kinematical analysis. 



knot ID 


i(°) 


v,o, (km s"') 


HH52A1 


58±4 


93±6 


HH52A3-A4 


57±2 


83±6 


HH52A 


62±4 


106±14 


HH52B 


63±3 


88±22 


HH52D2-D3 


67±4 


112±20 


HH52D4 


61±3 


83±11 


HH53 A 


81±1 


111±11 


HH53B 


85±1 


113±10 


HH53C 


83±1 


105+8 


HH53C1 


49±9 


76±21 


HH53E1 


57+3 


85±9 


HH53F2 


63±3 


100±15 


HH54 C2 


61±3 


51±11 


HH54 C3 


52±8 


24+9 


HH54G-G0 


55±8 


60-95±20 


HH54G1-G3 


53±8 


70±20 


HH54X1A 


66±6 


77±21 


HH54X1D-X4C 


64±7 


65±24 


HH54X3 


67±10 


50±15 


HH54X4A-X4B 


68±4 


124±23 



